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and Recovery of Submarine Cable 
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This paper is an attempt to jormulate a comprehensive theory vnth which 
the forces and motions of a submarine cable can be determined in typical 
laying and recovery situations. In addition to the fundamental case of a 
cable being laid or recovered with a ship sailing on a perfectly calm sea over 
a horizontal bottom, the effects of ship motion, varying bottom depth, ocean 
cross currents, and the problem of cable laying control are considered. Most 
of the results reduce to simple formulas and graphs. Their application is 
illustrated by examples. 
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Glossary of Symbols 

A Amplitude of harmonic ship motion 

Ci , C2 ,C2 Longitudinal wave velocity; trans- 

verse wave velocity in air and water 

(7q ^ Cj Transverse and tangential drag co- 

efficients 

d Cable diameter, also distance be- 

hind the ship at which the cable 
enters the lower stratiun 

D^ Dj, Noi-mal and tangential unit drag 

forces 

e Sidewise distance from the laid cable 

to the ship 

^A Extensile rigidity 

ji Ocean depth 

% ^ ^ Dimensionless ocean depth 

jj Hydrodynamic constant 

]j Inchned cable length from surface to 

bottom, also from ship to surface 

^^ Reynolds number 

p q Longitudinal and transverse devia- 

tional cable displacements 

p^ Qo Longitudinal and transverse ship 

displacements 
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Deviation from mean pay-out or 

haul-in rate 

Arc length and horizontal diHtance 
from the touchdo\vn pohit to the ship 

Dimensionless forms of S and X 

Time 

Dimensionless time 

Cable tension at an arbitrary point, 
at the ship, and at the bottom 

Dimensionless forms of T, 7\ and To 

Tp , Tq Cable tension due to longitudinal and 

transverse ship motion 

Ship speed, pay-out or haul-in rate 

Normal and tangential velocity of the 
water relative to the cable configura- 
tion 

Tangential velocity of the water rela- 
tive to a cable element 

Submerged and in-air unit cable 
weight 

a, ao Critical angle, approximate critical 

angle 

as Cable angle at the surface 

/3 Descent angle, cross current orienta- 

tion (Section 7.1) 



V, 


V, 


n 


■,Vr 


Vi 




w, 


Wa 



2 — sin a 
sin*' a 



Constant, also ascent angle 

e Slack 

Orientation of a cable element 

0^ ^ Spherical polar coordinates for the 

three-dunensional model 

K^ X Constants (see Appendix C) 
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Constant 

Constants 

Kinematic viscosity 

Rectangular coordinates for the three- 
dimensional model 

Mass densitj'' of water 

Mass per unit length of cable in air 
and water 

Deviation from the stationary angle, 

also angle between f axis and V 
(Section 7.1) 

I. INTRODUCTION 

In the summer of 1857, the first attempted laying of a transatlantic 
cable ended dismally when, after only a few hundred miles had been 
laid, the cable broke and fell into the sea. Although fouling of the pay- 
out gear caused by a negligent workman was the principal suspected 
reasons for the failure, its occurrence aroused great interest in the de- 
tailed dynamics and kinematics of the laying of submarme cable, and 
leading British scientists such as Kelvin and Airy pubhshed analyses of 
this problem in late 1857 and early 1858.^- ^' ^' '*' ^ 

However, after this initial activity, interest in submarine cable dy- 
namics and kinematics evidently waned for there appear only sporadic 
subsequent investigations in the literature.^' ^' ^' ^' '" Further, the re- 
sults of the early and subsequent analytical investigations have been, 
by and large, little utilized in cable laymg and recovery practice. One 
can conjecture several reasons for this. For one, because the early ana- 
lytical work was done before the advent of modern hydrodynamie theory, 
it did not rest on a secure base. Thus, as late as 1875, one finds vigorous 
debate over the nature of the tangential resistance of water to the cable.^ 
F'or another, the results of the analyses could not all be expressed in 
terms of elementary functions and required the numerical evaluation of 
some definite integrals. In the 1850's this was a tedious and laborious 
process. However, these are probably secondary reasons. For, after 
another failure in the early summer of 1858, a transatlantic cable was 
successfully laid in August of that year. The mechanical problem of de- 
positmg a cable was thus proved surmountable without complicated 
mathematical analyses, and the marriage of analysis and practice was 
never fully realized. 
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However, a present-day submerged-repeater transoceanic cable is a 
delicate and expensive transmission system. Reducing the amount of 
cable deposited by as little as one per cent can result in a substantial 
.saving in the first cost of such a system. Its repair is a costly operation 
requiring the sustenance of an ocean ship and its crew. Therefore, it is 
important to lay the cable without wasteful excess and with minimum 
chances for failure after laying. Further, it is important that repair, if 
necessary, be as efficient as possible. To accomphsh these things, an un- 
derstanding of the dynamics and kinematics of cable laying and recoveiy 
is essential. 

The purpose of this paper is to provide some of this understanding in 
as sti-aightforward a way as possible. To this end concepts and results 
are stressed in the main part of the paper, mathematical details being 
given in the appendices. Moreover, we hope to show that the results of 
the analysis can provide a numerical basis for decision making in many 
of the laying and recovery operations. Mo.st of these results can he ex- 
pressed in the foi-m of simple formulas and graphs. Several numerical 
examples are included to illustrate concretely how the results can be 
applied in practice. 

The general plan of the paper is to proceed from simple to more re- 
fined models of the laying and recovery processes. Thus, we discuss first 
what we have (tailed the two-dimensional stationary model. This model is 
appropriate for laying and recovery on or from a perfectly flat bottom 
while sailing on a perfectly still .sea. As a preliminary to this discussion, 
we consider in .some detail the hydrodynamic behavior of typical deep 
sea submarine cable. We then take up the effects of the ship motions 
which are induced by wave action and the effects of a bottom of varymg 
depth. These considerations are followed by a short discussion of the 
problem of controlling the cable pay-out properly during laying and the 
associated problem of the accuracy of the present taut wire method of 
determining ship speed. Finally, we consider the three-dimensional sta- 
tionary model and the effects of ocean cross currents. 

II. BASIC ASSUMPTIONS 

Our analyses, like most analy.ses of physical problems, are based on 
idealizations or mathematical models of the actual physical system. The 
extent of \'alidity of these models must be ultimately determined by ex- 
periment and experience. However we .-^hall try to give the reader an 
idea of when they are clearly applicable and when they are not. 

All (if the models we consider contain two basic idealizations, namely, 

(1) No bending stiffness in cable, i.e., it is a perfectly flexible string, 

(2) The average forward speed of the ship is constant. 
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Bending effects are caused by locally large curvatures, and are significant 
mainly where the cable leaves the pay-out sheaves and at the ocean 
bottom. However, for a cable with a steel strength member, bending 
even to the Hmall radius of the pay-out sheave typically does not ma- 
terially reduce the tension required to break the cable. Hence, in these 
cases we can expect an analysis based on the first idealization to give a 
reasonable idea of when cable rupture will occur. In laying, ship speeds 
are normally steady and, with the exception of the fluctuations caused 
by wave action which we consider later in the paper, the second idealiza- 
tion is reasonable also. In recovery, on the other hand, ship speeds are 
apt not to be steady, and the second idealization is more tenuous. But 
because of the very slow speeds usually employed, this idealization may 
in fact be meaningful in recovery as well. 

III. TWO-DIMENSIONAL STATIONARY MODEL 

3.1 General 

Assume that the cable ship is sailing at a constant horizontal velocity, 
that the cable pay-out or haul-in rate is constant, and that the drag of 
the water on the cable depends only on the relative velocity between 
the water and the cable. Further, assume that in a frame of reference^ 
translating with the ship the cable configuration is time-independent 
or stationary. This idealized model of the cable laying or recovery pro- 
cess we call the two-dimensional stationary model. 

This is the model which has been considered in the previous analytical 
studies.^'^" As the early investigators quickly pointed out, when the 
tension at the bottom of the cable is zero, the cable, according to this 
model, can lie in a straight line from ship to ocean bottom. During lay- 
ing, when slack is normally paid out, the zero tension condition actually 
occurs, and hence this case is of considerable practical miportance. 

The straight line can in fact be shown to be the only solution which 
can satisfy all the observed boundary conditions. This point is discussed 
in detail in Appendix A. That the straight line is a possible configuration 
can be seen from Fig. 1. In the vector diagram the velocity of the water 
with respect to the cable is resolved into a component Vn normal to the 
cable and a component Vt tangential to it. Associated mth Vn and Vt 
are normal and tangential water resistance or drag forces D^ and Dt ■ 
In the straight line configuration, the cable inclination is such that D^ 
just balan(;es the normal component of the cable weight forces. The 
situation is thus analogous to that of a chain sliding oji an inclined plane, 
with the forces Dn corresponding to the normal reaction forces of the 
plane. Summing forces in the normal direction, we get, therefore. 
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W COS a = D^■ , (1) 

while the HUinnialioii in the tangential direction gives for T^ , the ten- 
sion at the ship, 

Ts = wL sin a - DtL. (2) 

Here w is the weight per unit length of immersed cable, a is the cable's 
angle of incidence, Dn and Dt are the noi'mal and tangential drag forces 
per unit length respectively, and L is the inclined length of the cable. 
For most submarine cable used currently the force DtL is negligible 

and we ai'rive at 



w7j sin a — wh, 



(3) 



where h is the ocean depth at the cable touchdown point. Hence, during 
slack laying the cable tension at the ship is veiy nearly equal to the 
weight in water of a length of cable equal to the ocean depth. 







j^ r^°^' 




Fig. 1 ■ — Forces acting on a cable in normal laying. 

The straight-line .solution is the simplest and probably the most im- 
portant result to be obtained from the stationary two-dimensional model. 
We shall derive results for other important situations from this model 
also. As a preliminary, we study first, however, the nature of the normal 
and tangential di'ag forces D^f and Dt . 

3.2 Normal Drag Force and the Cable Angle a 

The resistance at sufficiently slow speeds to the flow of a fluid around 
an immersed body varies as the square of the fluid velocity. This relation- 
ship is usually written as* 



D^ = Cn 



pV/ d 

2 ' 



(4) 



* For towed stranded wire experimental verification of this relationship is 
reported in Reference 11. 
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Fig. 2 — Expcrimentiil imd thcoreticsil vjiriation of critical angle with towing 
velocity for cable No. 1. 

where Dn is the normal drag force per miit length, Co is the so-called 
drag coefficient, p is the mass density of the fluid, and d is the diameter 
of the cable. Vov the straight-line eonfigLU-ation, the vector diagram in 
Fig. 1 shows that 

V„ ^ V sin a. (5) 



Substitution of (5) and (4) into (1) yields m turn 

CnpV'd . 2 
w cos a = „ — sm a. 



(6) 



Equation (6) suggests how the value of the drag coefficient Co can be 
obtained experimentally. By towing a length of cable in water at a con- 
stant velocity, one can establish the straight-line configuration. The 
angle a can then be measured as a function of velocity, from which Co 
can be computed by (6). 

Figs. 2 and 3 show the results of such tests together with plots of (6) 
for the indicated values of Co - These results are taken from an analysis 
by A. G. Norcm of experimental data obtained by H. N. Upthegrove, 
J. J. Gilbert, and P. A. Yeisley. The properties of these cables are listed 
in Table I.* To eliminate end effects different lengths of cable were towed, 

* Cable No. 2 is very similar to present type D transatlantic telephone cable. 
For engineering calculations, type D can be con.^idercd the same as cable No. 2. 
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Table I — Properties of Cables No. 1 and No. 2 



Cable 


No. 1 


No. 2 


Wt. ill wiitcr (Ib.s/ft.) 


0.75 

0.243 
Polyethylene 
Smooth 


1,25 
705 


Outer covering 


Tu.r impregnated jute 
Rough 

4 X 10" lbs 

1.2 X 10" lbs 


EA (twist restrained) 

EA (Iwist iinrcslrained) 



as h iDfiicated by the plotted experimental points. It i.s seen that (6) 
gives a good fit to the experimental data over the eutire velocity range. 

If the cable has a smooth exterior, an estimate of the drag coefficient 
Ci, can be computed from published values of resistance to flow about 
an innncrsod cylinder, This computation is described m Appendix B, 
where wc have also tabulated computed values of Co . For the smooth 
cable No. 1, the value of Cd obtained from Appendix B is 1.00 which is 
in fair agreement with the experimentally determined value of 1.11. 

Although the drag coefficient Cd is a fundamental hydrodynamic para- 
meter, it is not the most convenient description of the effect of the nor- 
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Fig. 3 — Experimental and theoretical variation of critical angle with towing 
velu{;itv for cable No. 2. 
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mal component of water velocity. For small values of the ineidence 
angle a 

cos a ^ 1, 

sin a fii a, 

and (6) is approximately 

where ao is the approximate value of a. The quantity {2w/CDpd) ia a 
constant for a given tiable. It brings together all the cable parameters 
which influence the magnitude of the incidence angle a. If the angle a 
for a given speed is determined accurately, as can be done in a towing 
test or with a sextant during over-the-stern laying, this quantity is easily 
computed. Because of its importance, we shall call it the hydrodynamic 
constant of the cable and denote it by H, namely. 

By virtue of (7) and (8) we may write 

aoV = H. (9) 

The constant H rather than the drag coefficient Co will be used from this 
point on. 

When the approximate relationship (9) is not valid, a can be obtained 

by solving (0). This gives 



where V is in knots and H in radian-knots. In terms of ao we obtain in 
turn 



cos a ^ 'l/ ^ + 4 "^o^ ~ 2 "" 



(11) 



This relationship is shown in Fig. 4, whei'e the incidence angle a is plotted 
as a function of the approximate incidence angle ao . It is seen that for 
a < 20° the difference between ao and a is negligible. 

Physically a as given by (10) is the angle the cable assumes in the 
straight-line shape for the velocity T^. However, m addition, (10) shows 
that a can be thought of as a dimensionless parameter which embodies 
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Fig. 4 — Variation of a with ao . 

both the hydrodynamic properties of the cable and the ship speed. Thus, 
even when the configuration is not a straight hne, we shall find it con- 
venient to express results as a function of the single parameter a, rather 
than as a function of the two parameters H and V . For this reason, 
following Pode,"' '^ we call a the critical angle. 

3.3 Tangential Drag Force 

Over the range of velocities encountered in laying and recovery the 
drag coefficient Co in equation (6) is essentially con.stant. However, the 
corresponding coefficient for the skui friction force associated ^\ith Vt , 
the component of flow along the cable, is not constant. For the cable of 
smooth exterior (cable No. 1), the expression 

Dr = CflpVCTd, (12) 

with Cf = 0.055/(A'^ff)" '^, was found to give good agreement i\-ith the 
expei'imental data, as is shoAvn by Fig. 5. Here Dt is the skin friction or 
tangential drag force per unit length; F, is the relative velocity of the 
water with respect to a cable element gi^^en for -straight-line laying by 



Vt = Vc — V cos a, 



(13) 



where Vf is the cable pay-out rate; p is the mass density of water; and 
Nr is the Reynolds immber defined as A'^^ ^ Vt L/v, where y is the kine- 
matic viscosity of water. The data of Fig. 5 are for 100 foot lengths of 
cable towed in fresh water at a temperature of 60°F. 
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From (12) wo find 



Dr - 



0.055 V 

-^^ P 



0,14 



K/- 



L"-' 



d. 



(14) 



This expression indicates that Dt for smooth cable depends on the in- 
chned length of the cable as well as the relative tangential velocity Vi . 
The form of (13) suggests that the flow tangential to a smooth cable is 
similar to flow past a smooth plate. In such flow a tui-bulant boundary 
layer develops wliich grows in thictjiess with distance from the leading 
edge, resulting in a length dependence of the type shown by (14). Since 
Fig. 5 refers to 100 foot cable lengths, (13) is probably not accurate for 
the magnitudes of L occurring in deep-sea laying, and should be used 
only to obtain the order of magnitude of C/ . 
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Fig. 5 — Experimental values of the tangential drag force for cable No. 1 
compared with those obtained by equation (12). 

For the cable with conventional jute outer covering, (cable No. 2), 
it was found that 



T>r = 0.01 F/ ■" 



(15) 



fits the experimental data obtamed by towing test (Fig. 0). 'Whereas in 
(15) the constant 0.055 is dimensionless, the constant 0.01 in this equa- 
tion has the dimensions necessary to give Dt in units of pounds per foot 
when Vi is in feet per second. We note that for this cable Dt is inde- 
pendent of the length of the cable. 
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Fig. 6 — ■ Experimental vahies of the tangential drag force for cable No. 2 
compared with those obtained by equation (15). 

The ratio of Dr to the tangential component of the cable weight force 
is given by Dt/w sin a. Equations (14) and (15) indicate that even for 
small values of a of the order of twelve degrees, Dt/w sin a is of the order 
of 6 per cent for relative tangential velocities V, of 1.0 feet/see. In many- 
situations Vi will be less than this value, and we can neglect Dt compared 
to w sin a. As we shall see later, this approximation greatly simplifies 
the differential equations of the two-dimensional stationary model. 

Historically, the question of the variation of Dr with Vi is of some 
interest. In one of the early papers of 1858 Longridge and Brooks^ as- 
sumed a velocity squared dependence. In 1875, W. Siemens^ attacked 
this assumption stating that Dt actually varied luicarly i\ith Vt . There 
ensued a debate in which many bitter words but few experimental data 
were displayed. In view of our present knowledge that the skin friction 
force, even in the simplest case of flow past a smooth plate, is the result 
of complicated boundary layer phenomenon, the existence of this con- 
fusion is not surprising. 

3.4 Sinking Velocities and Their Relationship to Drag Forces 

The studies of submarine cal)le forces in 1857 and 1858 preceded 
modern fluid mechanics by many years. To characterize the hydrody- 
namic forces acting on cable the early investigators used sinHng or 
sealing velocities rather than the more recently conceived drag coeffici- 
ents. The transverse sinking velocity lu was defined as the terminal ve- 
locity attained by a straight, horizontal length of cable sinking in water. 
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Similarly, the longitudinal sinking velocity Vs was the terminal velocity 
of a cable length sinking with its axis constrained to be vertical. If for 
a given cable the drag forces are functions only of velocity, the pa- 
rameters w, Vs , and v, , together with the laws of variation of the drag 
forces with velocity, completely define the hydrodynamic behavior of 
the cable. Since sinldng velocities are still used in submarine cable tech- 
nology, it is of interest to relate them to the more modern drag coeffici- 
ent viewpoint. 

In the case of transverse or normal flow around the cable, the variation 
of Dx \vith the square of the relative transverse velocity gives (Vn/us)^ =^ 
Dn/w, since at a transverse velocity equal to the sinking velocity the 
unit transverse drag force is w. Substituting for Dn from (4) we find 

Thus, the transverse sinking velocity u^ is identical with the hydrody- 
namic consant H. We can therefore alternatively write the approximate 
relationship (9) as 

aoF- u,, (17) 

where ao is in radians and u, and V are in knots. 

For the tangential or skin friction flow along smooth cable, the sinking 
velocity concept is inadequate because the unit tangential drag force 
Dr varies mth length as well as with the relative tangential velocity 
Vt . However, for cable with the conventional jute exterior (cable No. 
2), we have (Vt/VsY'*^ = Dt/w and from (15) the vertical sinking velocity 
Vb is Vs = (46. Iw)^^^*^, where Vg is in knots. 

We note in passing that the cable does not, as is sometimes supposed, 
sink vertically to the bottom at the transverse sinking velocity u^ ■ 
Actually, the term "vertical cable sinking rate" is ambiguous. There are 
in fact two vertical sinking rates which may be important. Although 
both these rates are normally approximately equal to u^ neither is identi- 
cal to it. 

Relative to the earth, the resultant velocity Vr of a cable element has 
two components : a horizontal component of the magnitude of the ship 
velocity, and a component inclined at the angle a of the magnitude of 
the cable pay-out rate V^ . These are shown in Fig. 7. The component 
Fvert of Vr , glvou by Fyprt ^ Vc slu «, Is the rate at which a cable 
clement sinks vertically. For a laying depth h, the time r it takes for a 
cable element to sink to bottom is therefore r = h/Vc shi a. This time 
would, for example, tell one how long it takes a lightweight repeater, 
integral with the cable, to reach the ocean bottom. 
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Vt 




P' 
Fig. 7 — Illusti'iition of vortical cable siuking rates. 



On the other ha,iid, consider the intersection of the cable conjigtiration 
with a vertical line (Fig. 7). In the time t, as the ship sails a distance 
Vt, the intersection moves from A to ,1', a distance Vt tan a. Hence the 
cable configuration in this .sense sinks \'ertically at the rate V tan a, 
and the time d for the configuration to reach bottom in a depth h is 
6 = h/V tan a. One may be interested in how long it takes after the ship 
has passed over an ocean bottom anomaly P' (Fig. 7) for the cable con- 
figuration to reach the anomaly. This is just the time 6. 

Hence, the vertical sinking rates Vc sin a and V tan a can both be of 
interest. At the usual ship speeds, sin a ?t; tan a ;^ a ^ u,/V. Further 
Vc normally differs little from V. Hence, both these ratcw are indeed 
normally approximately equal to »« . 

3.5 General Solution of the Stalionary Two- Dimensional Model 

Assume that each cable element is traveling along the stationary cable 
configuration with the constant speed V^ . Starting at the ocean bottom 
let -s be the arc length along the stationary configuration. We define s to 
be positive in the direction opposite to the direction of travel of the cable 
elements. So, as Fig. 8 indicates, in laying, positive s is directed from the 
ocean bottom toward the ship, while in recovery the situation is reversed. 
We let 8 be the angle between the positive s direction and the direction 
of the ship velocity. 




^^m^^^^^^^;^^^^^:: 



Fig. 8 — Definition of coordinates for the two-dimensional stationiiry model. 
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Fig. 9 shows the forces acting on an element of the cable, with tension 
at the point s being denoted by T. The normal drag force per unit lengt-h 
Df/ may, by virtue of (4) and (5), be wTitten in the form 

It is necessary to introduce here | sin Q \ in order for Dn to have the proper 
sign for all B. We note, however, that if Y, ^ F, we have from (13) 

7, ^ V, - y cos e § 0. 

Hence in normal laying and recovery the unit tangential drag force 
Df is always in the positive s direction. 




wAs 
jri(,_ y — Diiigriim of forces acting on m cable elcn:ient. 

The forces acting on an element produce a centrifugal acceleration 
y/de/rfs. Hence, summing forces along the diret^tions i (tangential) and 
71 (normal), dividing by As and sending As to zero, we obtain 

(T - p^',^) - + ?^££^sin d I sine I - w cos d = 0, (a) 

^' ^ (18) 

— + Dr - w sin 6=0, (b) 

ds 

where p, is the mass density per unit length of cable. 

It is seen at the outset that e = a is a solution of (18a). It is in fact 
the important straight-hne solution which has been discussed in Sec- 
tion 3.1. 

If d 9^ a and Dt varies only mth Vt we may divide (18b) into (18a) 
and integrate to obtain the solution for T in the following form: 



(T - p.Vc') _ f jw s in g - Dt ) 

oa w(cos ^ — A sin 

CdP dV^ _ cos a 
2io sin- a 

where Tu is the tension corresponding to the angle ^o ■ 



, v^ Hr.-.) ^ f _ _(«Lsin_L_-^) ,ft (a) 

in - PcV,^) Joo w(cos ^ - A sin ? 1 sin ^ 1) ' 

^ Cop dV'^ ^ cos g 0^-j 

2io sin- a ' 
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At the cable touchdown point on the ocean bottom only two conditions 
are possible. If the angle 5 is not zero or w there, the cable tension T must 
be zero. Otherwise a finite tension would act on an infinitesimal length 
of cable, producing an infinite acceleration. Hence, either the tension T 
must be zero or the angle 6 must be zero or tt. The first case normally 
implies a straight-line configuration {see Appendix A), which has al- 
ready been discussed. In other cases, we define To as the tension at the 
touchdown point, and we let 6a be zero or tt, whichever is appropriate. 

If X, ij are coordinates in the translating (x, y) frame of a point along 
the cable configuration, then 

dx = ds cos 9, 

dy = ds sin 6. 

Combining these relations with (18a), we have 

{T - P.V.') 



jbo w{gos ^ — A sin f I sin f |) 

^ _ r' {T - P.O cos g 

iflo ■"'(cos ^ — A sin ^ I sin f |) 

_ r" {T - p.y/) sin e 



'Bo w{cos I — A sin f | sin ^ 1) 



d^, 


(a) 


dk, 


(b) (20) 


dl 


(c) 



Equations (19) and (20) are an integral representation of the complete 
solution of the basic two-dimensional model. In general, the integrals 
appearing in these equations cannot be evaluated m terms of elementary 
functions, and the solution must be obtained by numerical integration. 
For towing problems where the pay-out velocity is zero, Pode'^ has tabu- 
lated these numerical integrations using the approximation that Dt has 
certain constant values. However, in towing problems the direction of 
Dt is opposite to what it is in normal laying and recovery problems. Be- 
cause small magnitudes of i>r were used, these tables nevertheless usually 
give adequate results in laying and recovery situations as well. At the 
same time, for submarine cable problems, other approximations allow 
more convenient ways of evaluating the mtegrals of (19) and (20). 

For example, it is more accurate simply to assume that Dt is zero. As 
we indicated in Section 3.2, this approximation gives a negligible devia- 
tion from the exact solution if the relative tangential velocity Vi is small. 
Furthermore, in this situation we obtain from (18b) 

dT . ^ dy 

-r— = w Sin B — w ~r- , 
as ds 
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and hence the tension at the ship T, is very nearly 

T,= T,-\- wh, (21) 

where h is the depth at the touchdo-mi point. Thus, if the tangential drag 
force is negligible, the tension at the ship is essentially the bottom ten- 
sion plus wk, regardless of the nature of the normal drag forces. This is 
in fact a form of a well-known theorem which, as we shall see in Section 
7.1, applies in the three-dimensional case as well. 

In the next sections we make further simplifications of the general 
solution for the specific cases of laying and recovery. 

3.6 Approximate Solution for Cable Laying 

On long cable lays ship speeds are normally of the order of 4-8 knots, 
with accompanying values of the critical angle a of the order of 10°-30°. 
For these small values of a, the assumption of zero tangential drag to- 
gether with some mathematical approximations allow further simplifi- 
cations of the general solution. These simplifications are derived in de- 
tail in Appendix C; here we indicate the results. The angle 6 which the 
configuration makes with horizontal is closely given by 

, e , gf 1 - [n/(fo + y)r T 

tan 2 = tan - — — ^22) 

^ ^[l + [7'o/(ro + S)]nan4|J ' ^ ^ 

where y and To are dimensionless depth and bottom tension defined by 

y = y/K 

Ta = To/wh. 

Here we use the cable angle a in the sense of Section 3.2, namely, as a 
parameter characterizing the hydrodynamic cable properties and the 
ship speed. The constant y is in turn defined by 

(2 ~ sin a) 



y = — ^^ . (23) 

sm^ a 

For small a, tan* (a/2) is negligible and y is large. Further 

Hence, the denominator in (22) is very nearly unity and d approaches 
the critical angle a at small values of y, even for relatively large values 
of Ta of the order of three or four. Thus in the laying case, the cable 
configuration is very close to a straight line except for a short distance 
at the ocean bottom. 
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In Appendix C it is further sho\vn that for small a 

S = L-\- kTo/w, 

X = Lcosa -\- 'hTo/w. 



(24) 



Here S and X are the distance along the cable and the horizontal dis- 
tances respectively from the touchdowTi point to the ship (Fig. 11), 
L and L cos a are the corresponding distances for straight-line laying at 
the same ship speed, and k and \ are functions of the critical angle a 
which are plotted in Fig. 10. To illustrate the use of (24) we consider the 
following. 

Example: Cable No. 2 is being laid without slack onto a rough bottom 
from a ship moving at six knots. If the pay-out rate is decreased so the 
slack is 1 per cent negative, what is the subsequent rise of the tension 
with time at the ship? 

This is really a tran.sient problem. However, we shall try to get an idea 
of the average behavior of the cable by assuming it passes through a se- 
quence of stationary configurations. Also, we assume that because of 
the rough bottom there is no slippage of the cable along the ocean floor. 

If 5 is the amount of negative slack and V the ship speed, then in a 
time t an amount V{1 — S)t of cable will have been paid out. This amount 
plus the inclined length L will equal the amount contained in the curve 
AOC (Fig. 11). We then have 

L -j- V(l - 6)t= S -\- Yt - (X - L cos a). (25) 

Substituting (24) uito this equation and solving for To we find To = 

(w/(K — K))5Vt and that by (21) the tension at the ship is given by 



T, = wh-\- 



w 



\ ~- K 



5Vt. 
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Fig. 10 — Variation of k and X with the critical angle. 
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Fig. 11 — Cable geometry at a time t after the oiis^et of negative slack. 

For cable No. 2 a ship speed of sLx knots corresponds to a = 1 1 .7 de- 
grees. By Fig. 10, this corresponds to X - « - 1.4 X 10"^ Also w = 
0.705 lbs/ft by Table L Wc get therefore 



T, ^ wk-\- ;^ooo 



lbs 



nun. 



Thus, according to this calculation the tension in this example rises at 
the extremely rapid rate of 3000 Ibs/min. We note also that the rate of 
tension rise is here independent of the depth h. 

In the model which has been postulated, the cable is inextensible ; that 
is, it is assumed not to stretch under load. Because the difference between 
the lengths of AOC and the sum of the linear segments AD and DC (Fig. 
11) is small, one might suspect that the effect of cable extensibility in 
the present example is important. We can account for this effect in a 
crude way by assuming that the curve AOC has an additional length cor- 
responding to the stretching caused by the load To acting over the length 
L. For a cable made of a .single material, the stretching would be ToL/EA , 
where E is the Young's modulus of the material and A is the cross-sec- 
tional area of the cable. In analogy to this we denot^he extensile rigidity 
of the cable by EA, using the bar to indicate that EA is actually a single 
number obtained directly by measuring the extension of a length of cable 
loaded in tension. With this notation (25) becomes 



T,L 
L+ F(l - 5)t+YA 



= .S' + VI - {X - L cos a), 



and repeathig the previous computatioit we find 

w 



T, = wh + 



wh 



EA sin a 



+ X - K 



75(. 



It is to be noted that in this computation, unlike the iiiextensiblc case, 
the rate of tension rise depends on the depth h. 



DYNAMICS AND KINENL^TICS OF SUBMARINE CABLE 1149 

For (conventional holinally armored cable, one cannot define a single 
extensile rigidiiy hecauwo of coupling between pulling and twisting. Thus, 
how such a cable extends under tension depends on how it is restrained 
from twisting at the ship and at the ocean bottom. Instead of trying to 
determine these end restraints, we consider the limiting cases of no re- 
straint and complete restraint to twisting. Data supplied by P. Yeisley 
indicate the values of PJA for cable No. 2 in these conditions to be those 
given in Table I (Section 3.2). If we take h = 6,000 and 12,000 feet, we 
find with these values that 

h = 6000 feet: 

T, = wh + 220 (Ib/min); (twist mirest rained), 

= wh + 640 (lb/min)< (twist restrained), 

h = 12,000 feet: 

T, = wh -\- 120 (lb/min)i (twist unrestrained), 

= wh -\- 360 (lb/min)f (twist restrained). 

Comparing with the inextensible computation, we see that the extensi- 
bility markedly reduces the rate of tension build-up. Nevertheless, even 
for the case of no restraint to twisting at a depth of 12,000 feet the rise 
rate is a relatively high 120 Ib/min. Hence, at least over a I'ough bottom, 
the tension would quickly indicate the onset of negative slack, although 
the sensitivity of this indication would decrease with increasing depth. 

3.7 Approximate Solution for Cable Recovery 

Fig. 8 illustrates how cable is in present practice recovered from the 
ocean bottom. The cable is in front of the ship as it is brought in over the 
bow, and the ship pulls itself along the cable. In this process the cable 
tends to guide or lead the ship directly over its resting place on the ocean 
bottom. 

It is clear that during recovery by this procedure the tension at the 
ocean bottom is not zero and the cable configuration is not a straight line. 
Furthermore, in this situation the normal component of the water drag 
force Dn pushes down on the cable instead of buoying it up as in the case 
of laying. This in turn implies a higher tension at the ship during re- 
covery than during laying. 

If the tangential drag is neglected, the tension at the ship 1\ during 
recovery is given in dimensionless form by (.see Appendix C) 

^'^ -~ ^ r„. 2 cos a + cos a, ~|"^ ,.,.. 

= = tan a , (2G) 

is L 1 ~ <'0S a CO.S aj 



1150 THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1957 



wh 





































5 




w 


v 


V 






^ 


^^ 




















\v 


A 


s 


\, 










^^ 


^ 


^^ 








4 






\\ 


\ 


\. 


s 


\ 












«S=40'' 


"^ 


— 






\ 


.y 


N 


\ 




^ 


^^ 
















3 








N 


\^ 


\ 


\ 


^ 






"^ 


^ 


. 


50° 














\ 


\^ 


k 


^ 


^ 


-^ 










' 


— ■ 


2 












N 


\ 


-.^ 


"^ 


^ 




^ 


^ 


__60 = 


.. 


















^ 


■-- 


Zz 


c^ 


" 


■^ 


,^70° 







1 
























"^ 


-- 


3;^ 


-— 



30 40 60 60 

OC IN DEGREE5 



Fig. 12 — Viii-iiLtion of the tension factor for recovery with the critical anglea. 




2 3 4 5 

SHIP SPEED IN KNOTS 



Fig. 13 — Vivriation of exact and approximate tension factors for recovery 
of cubic No. 2. 
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where T, is the tension factor defined by T^ = Te/wh and y is given by 
(23). Equation (26) is plotted in Fig. 12 in the form of T, versus a for 
various surface incidence angles as (Fig. 8). It is seen that the recovery 
tensions are in fact considerably higher than the laying tension of approx- 
imately wh. 

To illustrate the smallness of the error of neglecting the tangential 
drag force in this computation, we have plotted the approximate and 
exact curves of 7^ versus ship velocity for cable No. 2 in Fig. 13. The 
dotted curves have been computed from (26), while the sohd curves have 
been obtained by substituting Dt from (15) into (19) of the general solu- 
tion and integrating numerically.* (The curve labeled Shea's recovery 
method is discussed in the next section.) 

The distance along the cable S and the horizontal distance .Y from the 
touchdown point to the ship cannot be expressed in a .'dimple form as in 
the case of laying. However, they can be obtained by numerical integra- 
tion from (20). The results of this computation for Dt = are shown 
in Figs. 14 and 15. 

How Figs. 12, 14 and 15 can be used is illustrated in the following ex- 
ample. 

* The stiindfird form of Simpson's rule was used for all the numerical integra- 
tions mentionod in the imi)or. In eiich case the interval of integration wan chonen 
fine enough to obtain at least three significant figure accuracy. 



\ 
































































y\ 


\ 


s 




























^ 


s\ 


V 


-^ 


^ 
























^ 


\^ 


x 

v 


^ 




^ 


—- 


— 


— 





_____ 




Kg = 40° 






^ 


^ 


— ^ 


^ 


^ 


■^ 


-- 


■— 


. 




■ — ■ 






50" 










^ 




^-■ 


^^ 


-^ 


■ — 


— 


. 




"= 


_60° 
70° 












"^ 


^^ 


::::; 


---- 


' ' 


— - 


— 


— 


. 


■— 


C 


















— 


— 


~- 


— 


-— 


so" 




























90" 


^ 


-^ 



40 50 60 

Oi IN DEGREES 



Fig. 14 — Variation of the horizontal distance to the touchdown point during 
recovery with the critical angle. 
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Fig. 15 — Variation of the distance along the cable to the touchdown point 
during recovery with the critical angle a. 

Example: A cable weighing 0.7 lb/ft in sea water and having a hydro- 
dynamic constant H of 70 degree-knots, is to be picked up from a depth 
of two thousand fathoms. If the ship speed is one knot what is the cable 
tension at the ship for surface angles a, of 40^ 60°, and 90"? How far in 
front of the ship and how far along the cable will the touchdown point 
be for these values of «« ? 

As indicated by (9), an H value of 70 degree-knots together with a ship 
velocity of one knot yields ao = 70 degrees. Fig. 4 yields in turn a = 
60 degrees. Entering Fig. 12 with this value of a, we can obtain TJwh. 
Ill this example the wk tension for a depth of two thousand fathoms is 
8,400 lb, and hence the values of T^/wh and T, are as follows: 



at 


T,/wh 


r. 


m° 


4.85 


40,700 lbs 


60° 


2.58 


21,700 lbs 


90° 


1.53 


12,900 lbs 


From Fig. 14 we get in 


turn for the horizontal distance from the ship to 


the touchdown point 






Olf 


X/h 


X 


40° 


2.65 


5300 fathoms 


60° 


1.5G 


3 120 fathoms 


90° 


0.66 


1320 fathoms 
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Finally, from Fig. 15 we get for the distance along the cable to the touch- 

doT\ai point 



as 


S/h 


8 


40° 


2.88 


5760 fathoms 


60° 


1.95 


3900 fathoms 


90° 


1.33 


2660 fathoms 



3.8 Shea's Alternative Recovery Procedure 

The high tensions which result in the usual recovery operation require 
slow ship speeds of the order of one knot or less if the cable is not to be 
broken. One wonders if it is possible to mitigate these tensions and thus 
speed the recovery process. J. F. Shea discovered that this can theo- 
retically be done by allowing a, to exceed 90°, thus establishing the 
straighthne configuration (Fig. 16). As in laying, the normal drag forces 
in this scheme support the cable, rather than push down on it as in con- 
ventional recovery. However, in contrast to the laying situation we have 
Vt = Vo + V COS a. Thus Vi is the sum of V^ and V cos a instead of 
their difference and Dt is not necessarily negligible. Furthermore, the 
direction of Dt is now such as to increase rather than decrease the tension 
over the wh value. Hence, instead of (2), a summation of forces along the 
cable yields Ta = wh -\- DjL, and the tension at the ship can be consid- 
erably higher than wh. A curve of 7", as a function of ship speed for cable 
No. 2 is shown in Fig. 13 with the label "Shea's recovery method". This 
has been computed for the case of haul-in speed equal to ship speed by 
means of the above equation and (15). It is seen that the tensions com- 
puted for this method of recovery, at least for the cable No. 2, are never- 
theless considerably .smaller than those which occur in the usual recov- 
ery procedure. It would .seem that the straight-Hne recovery technique 
could fruitfully bear further examination, especially for application to 
the recovery of long stretches of cable. 




^m^m^'^wM^m*^^^^^-^.^ 



Fig. 16 — The present and Shea recovery methods. 
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IV. EFFECTS OF SHIP MOTIONS 

4.1 Tensions Caused by Ship Motions 

In the basic stationary model a perfectly calm sea is postulated. 
However, in reality, wave action gives rise to a random motion of the 
ship which in turn induces variations in cable tensions around those cor- 
responding to the basic model. 

To analyze this effect, we assume that the mean forward velocity of 
the ship and the mean pay-out or haul-in rate are constant and that the 
mean tension at the ship and the mean direction of the cable as it enters 
water are those given by the stationary model. In a reference frame mov- 
ing with the mean velocity, we resolve the ship displacement into a 
longitudinal component Pa (Fig. 17) along the mean or stationary direc- 
tion and a transverse component Qo perpendicular to the stationary di- 
rection. 




Fig. 17 — Longitudinal and transverse components Po and Qo of the ship 
displacement. 

Intuitively, one might expect the tensions caused by the transverse 
displacement Qo to be negligible compared to those caused by the longi- 
tudinal displacement Po ■ An analysis we have carried through in fact 
yields this conclusion. Because of its complexity and length, this analysis 
and the model upon which it is based are given in Appendix D. The re- 
sults for the case of harmonic variation of Qo with time indicate, at least 
for cable No. 2, that the tension associated with the transverse com- 
ponent Qo is indeed negligible for all except ship motions so extreme as to 

rarely occur. 

In addition, this analysis indicates that for the transverse disturbance 
Qo , the amplitude of the responding transverse cable motion decreases 
exponentially after the cable enters the water because of the damping 
action of the water drag forces. The "half-hfe" distance for cable No. 2, 
that is, the distance along the cable at which the amplitude of a harmonic 
transverse motion is damped to one-half its surface value, is plotted in 
Fig. 18 as a function of the period of the motion for various depths h and 
ship velocities V. The striking feature of these figures is the rapidity of 
this damping. The analysis thus shows for cable No. 2 that the eflect of 
a transverse disturbance penetrates only a short distance into the water. 
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Fig. 18 — Variation of half-life distance of cable No. 2 with tho period of 
ship motion. 

As far as cable tensions are concerned, the important ship displacement 
then is the longitudinal component Po , directed along the stationary 
direction of the cable. The analysis of Appendix B leads to the basic 
one-dimensional wave equation 



b -p 1 ct p 



= 



for the description of the longitudinal motion. In this equation p is the 
deviation in longitudinal displacement from the mean pay-out or haul-in 
displacement, and the remaining symbols are defined as (Fig. 17) 

X = distance from the mean .ship position along the stationary cable 

configuration, 
t = time, 

ci^ = EA/pc ■ 



The additional tension Tp due to ship motion is in turn given by 



Tr 



^-RA^Jt 



dx 



(28) 



As in the example of Section 3.G, we have again assumed that by using 
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in (28) the limiting values of"^ obtained by complete restraint to twist- 
ing and no restraint to twisting during pulling, one can obtain bounds 
on the actual displacements and tensions. 

The solution of (27) under arbitrary boundary conditions can be ob- 
tained from standard textbooks. Probably it is most representative to 
assume the cable is semi-infinite. That is, although damping of the cable 
is normally so small that we neglect it in (27), we may assume, because 
of the cable's great length, that the damping is sufficient to cause com- 
plete decay of a disturbance initiated at the ship, and that such a dis- 
turbance is not reflected from the ocean bottom. Under this condition 
the additional tension Tp is given by 

T,= -VH^f, (29) 

where P = Pa ^- Pi "'ith dP^/dt being in turn the inciement in pay-out 
rate or decrement in haul-in rate from the mean. For cable No. 2, Table I 
(Section 3.2) indicates that 

^/ITApo ^ ^^^ lb/ft /sec (twist unrestrained), 

= 400 Ib/ft/sec (twist restrained). 

Two examples will make clear the application of (29). 

Example 1: Steady-State Laying or Recovery in a Regular Seaway. 

Assume that in a frame of reference traveling at the mean horizontal 
ship velocity ship surging (to and fro forward motion) is zero and the 
combined heave and pitch motion is normal to the ocean surface and is 
given by 

W = A sin27r-. 

T 

If the period t is seconds and the amplitude A is 15 feet find for cable 
No. 2, 

a) {Tp)^^,, for laymg at a constant pay-out rate and at a ship speed of 
6 knots, 

b) (7'p)„ax for recovery at a constant hau!-in rate and mth a surface 
incidence angle of 60°. 

In both cases (a) and (b), the deviation Pi in pay-out or haul-in rate is 
zero, hence P = Po = Tf sin a, , and (rfPM)nmx = {2-K/r) A sin a, . Since 
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nable No. 2 has a hydrodynamic constant H of 70 degree-knots, we have 
in case (a) 

a, = 11.7 degrees and ( -- ) = 2.12 ft/sec. 

\ di /max 

From (29) we get therefore 

(^p)mQx = 466 lbs (twist unrestrained), 
= 848 lbs (twist restrained). 
In case (b) we have «« = 60°, ((lP/dt)„tax = 13.6 ft/sec, and hence 
(Tp)max = 2,990 lbs (twist unrestrained), 
= 5,430 lbs (twist restrained). 

During recovery by conventional methods the surface incidence angle 
a, is in general much larger than that which occurs during laying. The 
above example points up that one can expect correspondingly larger 
ship motion tensions during recovery than during laying in the same 
sort of seas. Since the stationary tensions are also much larger during 
recovery, recovery is the condition for which the strength of the cable 
should be designed. 

In this example we have considered a regular seaway, something which 
does not exist in nature. Recent work in the application of the theory of 
stochastic processes to the study of ocean waves and ship dynamics 
promises to develop into a realistic description of the behavior of ships 
at sca.^^ When such a description becomes available, we shall be able to 
obtain a better estimate of the magnitudes of .ship motion ten.sions. 

As far as data presently available are concerned, the maximum storm 
condition vertical velocity at the bow or stern recorded by the U.S.S. 
San Francisco during lier research voyage of 1934 was 22 feet/sec. 
Since this ship was roughly the size of a cable ship such as the H.M.S. 
Monarch, this figure might indicate the order of the maximum velocities 
to be expected in cable practice. In terms of our example, for .six knot 
laying this vertical velocity would imply 

Tp = 980 lbs (twist unrestrained), 

= 1 ,780 lbs (twist restrained) . 

For recovery at a surface incidence angle of 60°, it would imply in turn 

Tp = 4,200 lb (twist unrestrained), 

= 7,600 lb (twist restrained). 
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However, it is to be cautioned that these nimibers are merely indicative 
and might differ considerably from those which occur on a particular 
cable ship. 

Example 2: Brake Seizure 

While laying cable No. 2 at six knots in a perfectly calm sea, a sudden 
seizure of the brake occurs. What is the resulting initial rise in tension? 
Because of the cahn sea we have Po = 0. Therefore 

—r- = Pi= —V COS a, . 
at 

With the value of F ^ 6 knots and a corresponding a. of 11.7° (see Ex- 
ample 1) we have dP/dt = 9.9 ft/sec and hence from (29) 

Tp = 2180 lbs (twist unrestrained), 

= 3970 lbs (twist restrained). 

These values of Tp pertain only to the transient values occurring while 
the tension wave is being transmitted to the ocean bottom. If the seizure 
in this case occurred at a depth of three thousand fathoms, the time of 
transit to the ocean bottom would be only of the order of nine seconds. 
After reaching bottom our initial assumption of no reflection from the 
bottom would be violated and (29) would no longer hold. In reality the 
cable tension would continually increase at the ship and reversmg ship 
engines or some other action would be required to avoid rupture of the 
cable. 

V. DEVIATIONS FROM A HORIZONTAL BOTTOM 

5.1 Kinematics of Laying Over a Bottom of Vanjing Depth 

Ocean bottom topography is not everywhere flat and horizontal as 
postulated in the basic model. In the Mid-Atlantic ridge, for example, 
there exist bottom slopes of thirty or forty degrees. In other places sub- 
marine canyons with almost vertical sides have been found. Further- 
more, where the bottom is steepest it is most likely to be rocky and 
craggy since erosion tends to smooth out a sandy or muddy bottom. 
Therefore, it is important to know how the cable should be paid out to 
cover a bottom of varying depth. To help determine this, we extend 
here the stationary model to the case of a non-horizontal bottom. 

In Section 3.1 we indicated that if the cable tension at the touchdown 
point is zero the configuration according to the basic model is a straight 
line, regardless of how the cable is paid out. If the cable is paid out with 
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slack with respect to the bottom, the zero touchdown tension condition 
is fulfilled. ITence, under the proper slack pay-out, the cable geometry 
and, as we shall see, the cable kinematics are particularly simple. 

Essentially, we must consider two deviations from the horizontal 
bottom, namely, downhill or descent laymg and uphill or ascent laying. 
We consider these situations in turn, confining ourselves to bottoms of 
constant slope since any bottom contour can be approximated by 
straight-line segments. 

To cover a descending bottom, the cable pay-out rate must exceed the 
ship speed, Fig. 19(a). To cover an ascending bottom, the angle of in- 
cidence a of the cable, which as we have seen in Sections 3.1 and 3.2 de- 
pends only on the ship speed, must exceed the ascent angle 7, Fig. 19(b). 
Otherwise, the situation shown in Fig. 19(c) develops. Hence the critical 
parameters are pay-out speed and ship speed. 




(b) ASCENT (a > 7) 





(C) ASCENT (;'>Ci£) 

Fig. 19 — Cable geometry during straight-line descent and iiseent laying. 

During descent laying we see from Fig. 20 that in a time t an amount 
of cable equal to a -1- b must be paid out. Hence the required pay-out 
rate Vc is (a + 6)/i. But by straightforward trigonometry 



V, = 



a -\- b sin a + sin /? 



V, 



(30) 



t sin (a -\- 0) 

where ^ is the angle of descent and a is the straight-line incidence angle. 
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Fig. 20 — ■ Kinematics of straight-line descent laying. 
In accordance with usual terminology we define the slack e as 
€ = (T^. - V)/V. 



(31) 



We shall think of the slack as being composed of two parts: a Jill f, which 
is the amount of slack required for the cable to cover the bottom, and an 
excess, equal to e — /, which will normally be laid to provide a margin of 
safety. Substituting Ve from (31) into (30) we get the expression for the 
fill /. The result can be transformed to the form 



/ 






(32) 



The quantities a, /3 and / are normally all small quantities and we may 
make the approximations 



tan^ 



a 

2' 



t.n^.?, 



/ 



2+/ 



For a, < ^0° and / < 0.06, the error in each of these approximations 
is less than 3 per cent. Hence, with good accuracy we write (32) as 



f = ^ 
J 2 • 



(33) 



where a and /3 are expressed in radians. Further, we have from (9) that 
a in radians is very nearly a ^ H/V, where H is in radian knots. Sub- 
stituting this expression into (33), we get for the fill 
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Finally, using this expression for /in (31), we arrive at* 

V. - F = f . (30) 

Thus, the increment in required pay-out rate is essentially a function 
only of the descent angle /3 and is inde-pendent of ship speed. 

In the case of an ascending bottom for which a > 7, Fig. 19(b), posi- 
tive bottom slack may be obtained ^Wth a pay-out of less than the ship 
speed. The allowable decrement in pay-out rate is given by 

V -V. = ^, (36) 

that is, the same as the required increment for ascent laying. Likewise, 
the fill / in this case is simply f — ~ {Hy/2V). 

The only way to avoid the situation shown in Fig. 19(c) where a < y 
is to sail slowly enough to maintain an incidence angle a greater than 
the angle of rise 7. By (9), we have for most laying speeds aV ^ H. 
With good accuracy the condition a > 7 thus implies 

V <-. (37) 

7 

Therefore, for a given rise 7 the limiting ship speed is simply H/y. 

5.2 TiTne-Wise Variation of the Mean Tension in Laying Over a Bottom 
of Varying Depth 

In the (lases where the cable is paid out with excess onto a bottom of 
constant slope, the variation of the mean tension at the ship with time 
is easily computed. During descent laying the increase in depth 6 after 
a time t is by elementary trigonometry (Fig. 20) 

_ sin a sin /3 ,, 
^ ' sm (« + /3) ^- 

Hence, the rate of rise of the mean shipboard tension is 

dT wS sin a sin /3 ,, ,„„. 

— - = — = - — -, ; — r wV. (38) 

dt t sm (a + j3) 

Similarly, during an ascent lay for which the bottom is less steeply 



• Note that in (33), (34) and (35), //maybe replaced by the numerically identi- 
cal transverse settling velocity it, (see Section 3.4). 
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inclined than the cable (a > 7), the rate of decrease in shipboard ten- 
sion is 

dT sin a sin 7 -.r 

■ = —-. ^—- r WV. 

at sin {a — 7) 

Like negative slack laying on a flat bottom, the variation of tension 
with time depends greatly on the frictional characteristics of the bottom 
in cases other than the above. We therefore hmit ourselves to situations 
where the cable does not move with respect to the ocean floor. This case 
might be approximated by rough bottoms, where the cable might wedge 
itself between rocks. 

A nomograph gi\'ing a rough estimate of the rise of mean tension \vith 
time when a cable becomes completely suspended Is worked out in Ap- 
pendix E. In deri^dng this nomograph it is assumed that the cable takes 
on a sequence of stationary configurations. This assumption is probably 
reasonable if the time span of the tension rise is large compared to the 
time of passage of a tension wave from the ship to ocean floor and return, 
which as mentioned in Section 4.1 is of the order of 18 seconds. However, 
because of this assumption and others mentioned in Appendix E, we 
regard the tension variation computed by the nomograph only as a crude 
approximation. 

Fig. 21 shows the mean ship-board tension versus time computed by 
means of the nomograph for various slacks e, where e is defined by (31). 
The values which were used for the other parameters entering the calcu- 
lation were 

-^ = 3.1 X lo-^ 

EA 
a = 12°. 

Also shown on this curve is the tension rise computed for the case of 
laying down a vertical slope without excess. The rise for this ease is given 
by (38) with /3 = 90°. It is seen that as the slack e is increased the curves 
for a complete suspension approach the /3 ^ 90° curve. Indeed, it can be 
shown that under the assumptions made in computing Fig. 21 the ^ = 90" 
curve gives a lower bound on the tension rise Avith time in the case of a 
complete suspension. A tension rise rate greater than the /3 = 90° rate 
is thus an indication of unsatisfactory covering of the bottom. 

In the case of too rapid a ship speed resulting in a < 7 (Fig. 19c) re- 
straint of movement of the cable along a rough bottom ^\ould cause the 
tension on the high side of the crest to be zero. There would thus be a 
sudden drop in tension corresponding to the sudden decrease iu depth 
at the touchdown point after the cable was laid over the crest of the hill. 
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Fig. 21 — VariiitioL of tensiou with time when cable No. 2 is completely sus- 
pended. 

On the other hand, in the case of a frictionless bottom, the removal of 
the supporting water drag forces would cause the cable to seek a catenary 
equilibrium position on the low side of the crest. But in doing this, the 
cable would drag itself over the crest, with an accompanying increase in 
shipboard tension. 

Thus for the case of a bottom rise steeper than the cable inclination 
(a < 7) either an increase or a decrease of tension with time is possible, 
depending on the nature of the bottom. 



5.3 Residual Suspensions 

If the cable is not paid out rapidly enough, or if the ship speed is ex- 
cessive, the cable will be left with residual suspensions after it has been 
laid. To get an idea of the possible magnitudes of the tensions accom- 
panying these suspensions, we consider here some numerical examples 
pertaining to cable No. 2. As before, we assume for definitoness the ex- 
treme case of a bottom rough enough to prevent movement of the cable. 

In Fig. 22 is shown the profile of a 35 fathom (210 feet) increase in 
depth with a maximum slope of 45°. This profile was obtained from 
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fathometer records of the Mid-Atlantic Ridge provided by Professor 
Bruce C. Heezen of the Lament Geological Observatory. Laying domi 
this slope at a ship speed of six knots requires a slack of 8.5 per cent {see 
Section 5.1). If the slack were only 5 per cent, the successive cable con- 
figurations as calculated by the methods of Appendix E would be those 
shown in Fig. 22(a).* The cable would touch bottom after 2.6 minutes, 

t IN minutes: 






(a) DESCENT 



'T = 2710 LBS 




-990' *] 

1320' »| 

(b) ASCENT 



Tig. 22 — SuccesBive cable configurations during a 35-fathom descent and as- 
cent fay of cable No. 2 at a 6,000-ft depth witli an assumed 5 per cent slack and 
6-kuot ship speed. 

leaving a residual suspension with a half -span of 480 feet and a tension 
of 525 lbs. The mean tension at the ship would correspondingly increase 
by 525 lbs during the 2.6 minute time interval. In even moderately 
rough seas, this tension change could be obscured by the ship motion 
tensions. 

Consider this profile next to represent an ascending lay under a ship 
speed of six knots. Fig. 22(b) shows the initial {t = 0) and residual cable 
configuration. Because of the small incidence angle of the initial straight- 
line shape, the residual half-span of the catenary is a quarter of a mile 
(1320 feet) long, and the accompanying residual tension is 2,710 lbs, or 
roughly that which normally occurs in laying at a depth of f of a nautical 
mile. At the ship, there would be a decrease in the mean tension of 130 
lbs. corresponding to the 35 fathom decrease in depth. Again, a tension 
change of this magnitude would he difficult to discern because of ship 
motion tensions. 



* We have further taken the ratio wh /EA to be 3.1 X 10"" in tlii^ computation. 
However, the results are very insensitive to change in the wh /EA ratio. 
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If the above 85 fathom change occurred at a depth of say three thous- 
and fathoms, a very sensitive fathometer would be required to detect it. 
Thus, although complete restraint of cable movement along the bottom 
is an extreme and unlikely condition, the above examples indicate that 
long residual suspensions can occur with essentially no manifestation at 
the ship, especially in deep water. 

VI. CABLE LAYING CONTROL 

6.1 General 

We have seen that the mean cable tension at the ship reflects the 
amount of slack which is being paid out and how the cable is covering 
the bottom. However, in most cases this reflection is not sensitive. For 
example, the tangential drag force Dt varies with. Vi , the longitudinal 
velocity of the cable relative to the water. In theory, as (2) .shows, one 
can therefore determine the amount of slack being paid out from ship- 
board tension measurements. For cable No. 2, we have plotted in Fig. 
23 the variation of the mean tension at the ship as a function of slack for 
a ship speed of six knots and a depth of two thousand fathoms. At three 
per cent slack the tension is 8,240 pounds, while at six per cent slack it 
is 8,020 pounds, a difference of only 220 pounds. This amount of tension 
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Fig. 23 — Variation of Khipboiird tension with per cent slack for laying cable 
No, 2 at a ship speed of 6-knots in a depth of two thousand fathoms. 



could be easily obscured by the effect of ship motion. Thus, to measure 
slack a(;curately by relating it to cable tensions one would have to know 
the depth and cable parameters very precisely and, in addition, would 
need a very efficient filter to separate out the "noise" tension caused by 
ship motion. 

Similarly, it has been shown that re.sidual suspensions can occur with 
essentially no reflection in the tension readings at the ship. Hence, al- 
though tension readings can give a valuable check on how the cable is 
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covering the bottom, it would seem difficult for them to provide exact 
enough data for the control of cable laying. 

At the same time we have seen that if the bottom contour is known 
in advance, then for a given ship speed one can compute the required 
cable pay-out rate. Also with foreknowledge of the bottom, one can 
anticipate steep bottom ascents and decrease the ship speed accordingly. 
Such a purely kinematic attack on the cable laying problem would seem 
more fruitful than an attack which depends on measurements of ship- 
board cable tensions. 

Pos.'^ibly the simplest way of measuring the bottom contour is by 
means of a fathometer located at the ship. Since the cable ship is nor- 
mally far forward of the touchdown pouit of the cable, one could in 
theory obtain in this manner the required advance knowledge of the 
contour. In present practice, a taut piano wire is used to obtain the 
ground speed of the ship. We examine briefly the accuracy of this method 
in the next section. 

6.2 Accuracy of the Piano Wire Technique 

The taut wire is laid simultaneously with the cable, but under a con- 
stant mean shipboard tension. If the bottom is perfectly horizontal, the 
speed of the wire coincides with the ground speed of the ship. However, 
when the bottom depth is variable and the wire is laid up and down hill, 
the wire's pay-out speed deviates from the ship speed. By (31), it is seen 
that the error in the ship speed which is indicated by the wire is just 
equal to the slack e with which the wire is paid out. This slack, which 
may be positive or negative, can in turn be estimated by the methods of 
the previous sections. 

Consider the beginning (denoted by (1) in Fig. 24) and end (denoted 
by (2) in Fig. 24) of a doT^^lhill lay of the piano wive. As before we neglect 
the tangential drag force. Then, the condition that the tension at the ship 
remains constant gives, by (21), 

(To)i'hwh = iTo)2-\-wh, (39) 

where the subscripts 1 and 2 refer to the configurations at the beginning 
and end of the downhill lay. If e is the average slack or error of the piano 
wire during the descent, then (I + e}V is its average pay-out rate, and 
we have by Fig. 24, 

5, + (l+0Fi = ^i^ + 5,, (40) 

sm /3 
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where »S'i and S2 are lengths along the cable from the touchdown point 
to the ship. Also, from Fig. 2-4 

X, + VI = {]h - ^i)/tan )3 + A% . (41) 




Fig. 24 — Piano wire configurations at the beginning and end of a descent lay. 

Equations (39), (40), and (41), together with the general equations 
(Section 3.6) 



8 - 


h . To 
= -. — + K — , 
sin a w 


X = 


tan a w 



(24) 



allow one readily to solve for the average error i in the piano wire indi- 
cation of ship ground .'speed. The result is 



_ _ sin g + sin ;3 — K sin a s i n 13 _ 
' ~ sin (a + ^) - X sin a sin j3 



(42) 



For the small values of a and which normally occur during the laying 
of the piano wire, the terms \ sin a sin /3 and k sin a sin are neghgible. 
Hence, the average error e is thus very nearly 



sin a + sin j3 , 
€ = —. ; ■ ^s ^ 1, 



Sin {a + p) 

which, as (30) indicates, coincides with the amount of fill which would 
be required to lay downhill with the straight-line or zero touchdown 
tension configuration. Equation (43) is in turn closely approximated by 
(Section 5.1) 



^ ~ 2 2V' 



(44) 
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Similarly, for iiseent laying of the wire on a bottom which rises less 
steeply than the inclination of the wire (19b), we get 

sin a — sin 7 + K sin a sin 7 



« = 



sin {a — y) -\- \ sin a sin 7 



- 1. 



which is very nearly 



— ay 



Hy 
'2V' 



(4o) 



(40) 



Thus, in both the above cases, the error in the piano wire technique can 
be closely obtained by assuming that the configuration of the wire is a 
straight line during ascent and descent laying. This is not surprising 
since, as we saw in Section 3.6, the deviation from the straight-line con- 
figuration during piano wire laying is normally small. 

Because of its smooth exterior, the normal or transverse drag coefficient 
of the piano wire probably can be obtained from published curves for 
flow past a smooth right circular cylinder as shown in Appendix B. For 
typical 12 gauge (0.0290 inch diameter) piano wire, these curves yield 
a value of Co of 1.45 and an /f value of 25.0 degree-knots. However, these 
values of Co and H must be considered tentative until confirmed expe- 
rimentally. 

Knowing the wire's H value, we can compute the error of the ground 
speed caused by descent and ascent laying of the piano wire by means 
of (44) and (4G). The result of this computation for H = 25.0 degree- 
knots is sho^\Ti in Fig. 25. 

When the ascent angle of the bottom exceeds the incidence angle of 
the wire, suspensions result and the error cannot be computed without 




6 a 10 2 

SHIP SPEED IN KNOTS 



Fig. 25 — Error during descent and ascent laying of 12-gaiige piano wire. 
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knowledge of the frictional properties of the bottom. For an H value of 
25.0 degree-knots, (87) indicates that suspensions will occur for ship 
speeds V greater than V = 25.0/7, where V is in knots and the ascent 
angle 7 is in degrees. Hence, for a typical laying speed of (i knots, ascent 
angles greater than 4.2 degrees will cause suspensions of the piano wire. 
These magnitudes indicate that suspensions of the piano wire probably 
actuaUy develop in practice. 

It is seen from l''ig. 25 that for the usual small ascent or descent angles, 
the piano-wire technique is quite accurate, while for large bottom slopes 
it can be cou.siderably in error. Again, however, if the bottom contour is 
known in advance, these errors can be estimated in the cases plotted in 
Fig. 25 and therefore can be corrected for. In this manner, the piano 
wire could be improved to give accurate ground speeds in all two-dimen- 
sional situations, with the exception of the case of a suspension caused by 
a too .steeply ascending bottom. Such suspensions can be avoided only 
by maintaining a sufficiently slow ship speed. However, as seen by the 
small computed // value of 25.0 degree-knots, the ship speeds required 
to avoid piano wire suspensions on mieven bottoms are probably pro- 
hibitively slow. Hence, for steeply ascending bottoms it is likely that 
some other means of determining the ship ground speed is necessary. 

VII. THREE-DIMENSIONAL STATIONARY MODEL 

7.1 General 

Thus far we have assumed that the cable lies entirely in the plane 
formed by the ship's velocity vector and the gravity vector. Because of 
the symmetry of the cable cross-section, this assumption seems reason- 
able.* However in certain cases, as for example in the presence of ocean 
cross-currents, the assumption of a planar configuration is clearly un- 
tenable. We consider therefore the case where the cable configuration 
is not necessarily planar but is still time independent with respect to a 
reference frame translating with the constant velocity of the ship. In 
analogy with previous termmology, w^e call this the three-dimensional 
stationary model. 

Assume there is a constant velocity ocean current in each of a finite 
number of layers. Let the vector Vw denote the ocean-current velocity 
in a reference layer. In the stationary situation the velocity of the cable 

* Becaui^e of asymmetries caused by the helical armor wire or because of minor 

out-of-roiiiulnnsH, it is conceivable that a sidewise drag force might develop which 
would cause the cable to move out of the ship's velocity-gravity plane. For a re- 
port of experimental observations of such yawing in wire stranded cable.s, see Ref- 
erence 11. 
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configuration is everywhere the velocity of the ship, which we denote by 
the vector V. Hence the velocity V' of the water Avith respect to the 
cable configuration in the reference layer is 

r = V^-\- (-V) = 7. - F. 

Further, in this layer we choose a set of coordinate axes $, tj, f translating 
at the velocity V as follows: The f axis has the direction of —V', w^hile 
1} is measured vertically upward, and f is perpendicular to ?; and ^ so 
that the axes ^, rj, f form a right-handed system. A plan view of this 
configuration is shown in Fig. 26. We have denoted the angle between V 
and Vw by /3, while the angle between the ^ axis and V is denoted by ip. 
(The distances d and e refer to a subsequent section.) To describe the 
cable configuration with respect to the ^, rj, f axes, we use the spherical 
polar coordinates 6 and ^ show^l in Fig. 27. {Thet,u,~v vectors are dis- 
cussed in Appendix F.) 




Fig. 26 — Plan view of the coordinate system for the three-dimensional sta- 
tionary model. 

As in the two-dimensional case, we resolve the velocity of the water 
with respect to a cable element in the reference layer into a component 
Vff normal to the cable and a component Vt tangential to the cable, and 
associate with Vk and Vt the drag forces Ow and Dt . The resulting 
differential equations, which are derived in detail in Appendix F, are the 
following : 



-|- wA'(coh^ }p siii^ 6 -\- sin^ ^)^'^ cos ^|/ sin 6 — w cos 9 = 0, (a) 
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as 



(47) 



+ «;A'(cos" ^ sin 9 + sin^ ^py''^ sin ^p = 0, (b) 



— \- Dt — w sin 8 ^ 0, 

as 



i^) 



where A' = C„pdV''^/2, and V is the magnitude of f'. 

In addition, connecting the coordinates ^(s), 7j(s), and f(s) of a point 
s along the cable with the angles 6 and ^ we have the geometric relation- 
ships 

d^is) 



ds 

dyjs) 

ds 

rff(s) 
ds 



= cos cos ^, 



(a) 
= sin fl, (b) 

= —cos 6 sin ^. (c) 



(48) 



Two important general results follow from (47) and (48). For one, if 
the tangential drag force Dt is neghgibly small, (48b) substituted into 
equation (47c) yields upon integration 

7' = To + wv, (49) 

where To is the tension at ?; = 0. Hence, if t; is measured from the ocean 



t DIRECTION OF 
CABLE 




Fig. 27 — Definition of the spherical poliir coordiiiaLcs und iP iinil Iho unit 
vectors /, u, and v. 
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surface, and if at the bottom (77= —h) the tension is zero, the tension 
at the ship is essentially wh, regardless of the nature of the normal 
drag forces. Since in most laying situations for present cables, the tan- 
gential drag force can be reasonably neglected, this fact provides a con- 
venient over-all check on the laying process. That is, if the cable is being 
laid with excess, the tension at the ship for any stationary cable con- 
figuration, planar or non-planar, should be essentially wh. Any marked 
mcrease of tension over the wh value necessarily means the bottom ten- 
sion is non-zero and insufficient cable is beuig paid out. 

The second important result is derived in Appendix F. This result is 
that if the bottom ocean layer in our model is devoid of cross currents, 
and if the bottom tension is zero, then, for the boundary conditions which 
are normally observed, the cable configuration in the bottom layer is a 
straight Una. Further, this straight fine is in the plane formed by the 
ship's velocity vector y and the gravity vector. Hence, for example, in 
laying with excess in a sea which contains surface currents, the cable 
configuration in the lower, current-free portion will be a straight line in 
a vertical plane parallel to the resultant velocity of the ship. The laid 
cable will be parallel to the ship's path, but displaced a certain distance 
from it. Thus, because the lower portion is a straight line, our previous 
results about the kinematics of straight-line laying still apply. Only they 
now are pertinent to the displaced bottom contour rather than to the 
contour which lies directly beneath the ship. 

7.2 Perturbation Solution for a Uniform Cross-Current 

Cross-currents are commonly confined to a region near the ocean sur- 
face. It is of interest therefore to determine for such surface currents the 
distance e (Fig. 26) which the laid cable will be displaced from the path 
of the ship. In Appendix F we consider the problem for a cross-current of 
uniform but comparatively small velocity. In addition, we determine the 
distance d (Fig. 20) back of the ship at which the cable leaves the upper, 
cross-current stratum and assumes the straight-line configuration it has 
in the lower stratum. Let us assume for the sake of I'eference that the 
resultant ship velocity V is due east, and that the cross-current Vu is 
inclined at an angle ^ to the north (Fig. 26). The resultant velocity V 
of the water with respect to the cable in the surface stratum has the 
magnitude therefore of 

V = [(V - 7. cos ^f + (K. sin 0f]'', (50) 

and is inclined at the angle tp from due west, where 
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tan ip = Y^ :^ — - — z. . 

F - y,„ cos )3 



(51) 



Associated with V we have a critical angle a' which is given approxi- 
mately by H/V or exactly by (10) or (11). 

In terms of fp, a', and a the analysis of Appendix F yields the following 
values of d and e. 



= h ^tna I 



1 h -k' , 
1 _ __^ tan 



■■[ 



1 - 1 - 



h' 



r\ ctn-'d 



'D- 



(52) 



(b) 



where Aa = a — a is the difference of lower and upper stratum critical 
angles, h' is the depth of the upper, cross-current stratum and h is the 
total depth.* Curves from which d and e may be evaluated are given in 
Figs. 28 and 29 in dimensionless form. To illustrate their application we 
consider the following. 



- -6. 




Fig. 28 — Dist-iince the laid cable is offset from the .ship's paMi. 



* In Appendix F the equation of the .space curve formed by the cable in the 
upper stratum is also given. 
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Fig. 29 — Distance behind tiie ship at which the cable enters the lower stra- 
tum. 

Example. A ship is laying cable No. 2 at a depth of 6,000 ft at a re- 
sultant ground Kpcod of (i knots due east. There is a one knot cross-cur- 
rent GOO feet deep running 80° east of north. Fmd e and d. 

Here ^ = 60° and we obtain from (50) by a simple computation 
V = 5.57 knots. Also since for cable No. 2 H = 70 degree-knots, 

k'/h = 600/6000 - 0.1, 

a' = 70/5.57 - 12.3 degrees. 

By interpolation, we find from Figs. 28 and 29 



\tan a' h'J 



.tan a' h' / Aa 



Equation (51) yields in turn ^ — 0.156 radians, and Aa is given by 

70 70 
Aa = =^ — -^, = —0.6 degrees ^ —0.010 radians. 

Hence, wo have 

e = 3.6 h'<p - 2.7 X GOO X 0.156 - 253 ft, 

d = h'\ ' -. - 4.7Aa - 6000 [4.59 + 4.7 X 0.010] = 2800 ft. 

Ltan a J 
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APPENDIX A 

Discussion of the Two- Dimensional Stationary Configuration for Zero 
Bottom Tension 

We assume again that the tangential drag Dt depends only on the 
relative tangential velocity Vt , and we consider ma T, 6 plane the solu- 
tion trajectories of equations (18). These trajectories satisfy the equation 

dT _ (sin e - Dt/w) 
(id ~ 



{T - P.VI), 



(53) 



cos 8 — A sin 6 \ sm 6 

and arc periodic in 6 \vith a period of 27r. In Fig. 30 we have plotted the 
solution trajectories qualitatively for (a — tt) ^ 6 ^ (a + tt). It is seen 
that the trajectories are either the vertical straight hnes 6 = a,6 = azizT 
or they lie completely within one of four regions, labelled I, II, III, or 
IV, which are bounded by these vertical lines and the horizontal hne 
T = pcVc^. The trajectory 6 = a corresponds to the straight-hue lajdng 
configuration, while the trajectories 6 = a ± t correspond to Shea's 
straight-line recovery method. 

Examine now the trajectories in Regions II and III at a point of 
which r = 0. As J. F. Shea has pointed out, these trajectories all lie 
below the line T = PcVc. On the other hand, the trajectory 8 = a con- 
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Fig. 30 — Qualitative representation of the solution trajectories of the two- 
dimensional stationary model. 
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tains all values of T. Hence, according to the stationary model, the only- 
cable configuration for laying which has the value T — and values of 
T ^ pcVc is the straight line inclined at the critical angle a. The mag- 
nitude of pcVc is small. For example, for cable No. 2 paid out at six 
knots PcVe is roughly six pounds, and for conditions approximating 
stationary laying the observed tensions at the ship are in practice always 
many times the pcVc value. For such magnitudes of shipboard tensions 
and a zero bottom tension, the two-dimensional stationary model thus 
yields the straight line as the only possible cable configuration. 

However, the empirical fact that T > pcVc does not guarantee that 
the shipboard tension must be greater than pcVc- We might somehow 
contrive to lay at a zero bottom tension with T < pcV^ and with the 
cable in one of the non-straight line configurations of Regions II or III. 

Consider the cable configuration lying in Region II. From Fig. 7 it can 
be seen that the vertical velocity of a cable element is given by 

-3 = -"^vert = -K sin0, 
at 

where y is measured upward. Hence, of the possible trajectories for which 
the bottom tension is zero oidy those for which the bottom cable angle 
do is between zero and x correspond to cable laying. For region II, there- 
fore we need consider only the trajectories in the range ^ 5o < a at 
2*0 = 0. From (20c) the maximum value of y^ for these trajectories is 
given by 



Vm = 



f 



p^v: r I sin ^ 



w Jo [{cos f — A sin^ |) 



><^^PL"i£ w-CcosT,- Asin^T,)^^]/^^- 

Let (-Dr)m be the maximum value of Z>t , ^ tj < a. With Dt set equal 
to (Or)™ , the right-hand side of (54) gives an upper bound on ?/,„ . This 
substitution further allows one to evaluate the right-hand side of this 
equation in terms of standard integrals. The result yields the following 
upper bound on y^ '■ 

P.V.' 1 



y^ < 2.1 
where 



r = 



w 1 — r' 

w sin a ' 
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In general, this upper bound vdW be much less than the laying depth. 
For example, for cable No. 2 being laid with 6 per cent slack at 6 knots 
y,n < 12.5 feet. That is, the cable configurations corresponding to Region 
II do not reach the ocean surface. Hence these solutions of the stationary 
model do not in general satisfy all the required boundary conditions and 
can be discarded. 

Similarly, in Region III, the laying trajectories for which To = are 
in the range a < Bu ^ it. Consider those for which 6q < 7r/2. We get for 
these trajectories 






PcVi r J sin ^ 



w Jb, UAsin^f - cosf) 

(55) 



X exp 



L Jflo ^ (A sm^ V — cos t;) JJ 



where ijr/2 is the value oi y at 6 = ■ir/2. Let m be the minimum value of 
sin J7 — {DT{ri)/w) in the range a < 77 ^ 7r/2. If, as in the usual case, 
m is positive, we can obtain an upper bound on 7/^/2 by replacing 
sin Tj — (Dr(7i)/M) by m in the right-hand side of (55). By this means 
we find that 

^ P.y/ 2(1 + cos" a) 
w m tan a/2 

For cable No. 2 being laid with 6 per cent slack at G knots this relation 
yields y^/2 < 1,100 feet. So in the usual laying depths, which are many 
times greater than ?/,/2 , the configurations in Region III for which 
To = correspond to a value of 6 at the surface greater than t/2, or to 
cable being paid out in front of the ship during laying. It is doubtful 
whether such configurations would be stable and, at any rate, doubtful 
whether cable would ever be laid in such a manner. Hence, we conclude 
that these To = solutions of Regions II and III mil in general be 
mathematical curiosities, and that the only realistic laying solution of 
the stationary model for which the bottom tension is zero is the straight 
fine 9 = a. 

APPENDIX B 

Compulation of the Transverse Drag Coefficient and the Hydrodynamic 
Constant of a Smooth Cable from Published Data 

From (6) of Section 3.2 we obtain the relationship 



1178 THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1957 

-^ = (^Y'^ = „, (56) 

V COS a \L Dpd ' 

where H is the hydrodynamic constant. Also, we define the Reynolds 
number for flow normal to the cable in the usual way ; 

^^^dysma_ (57) 

V 

V being the kinematic viscosity of water. For smooth cable we now as- 
sume that ihe drag coefficient Ci> depends on Nr in the same way as for 
flow around a smooth cyhnder of infinite length. Experimental data for 
this relationship, namely, 

Co = CniNjd (58) 

are available in the Uterature and have been collated by Eisner. 

For a given velocity V, (56) through (58) represent three equations 
for the unknowns a, Co , and Nr . In general, the solution of these equa- 
tions depends on V, thus contradicting the assumption made in Section 
3.2 that Cp , and therefore H, are constants independent of V. However, 
for sufficiently large V we can expect the resulting a to be small so that 
Vcos a f^ 1. In this case (56) and (57) combine to give 

C. = ?^^„ (59) 

which together with (58) yields two equations for Co and Nr that are 
indeed independent of V. In laying, V \vill normally be large enough for 
this approximation to hold. 

Equations (58) and (59) are in turn easily solved graphically by find- 
ing the intersection on log-log paper of the curves, Co versus Nr , 
that these equations represent. Since p and v are properties only of the 
water, we see that Co is a function only of the product wd of the unit 
weight of the cable times its diameter. In Fig. 31 we have plotted the 
resulting values of Co for wd ranguig from 10"'' to 10 pounds. For this 
computation we have assumed sea water at 32°F with an assumed density 
of 1.994 slugs per cubic foot and a kinematic viscosity of 2.006 X 10~^ 
ftVsec. 

For other than large values of V, (56) through (58) can be readily 
solved if one interchanges the roles of a and V, that is, if one considers a 
as given and V as unknown. Equations (56) and (57) can again be com- 
bined in this case to give 
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Co = 



2wd cos a 1 



2 ' 



(60) 



and with wd cos a a. known number, (60) and (58) can be solved for Co 
and Nn as before. Thus, one can obtain Cd from Figure 31 by merely 
reading wd cos a rather than wd on the abscissa. Knowing Cd one can 
solve for V from (56) . 



3.0 



2.0 



wd IN POUNDS ( ) 

2 4 6 O'O"'^ 2 4 O'O-* 2 46 8 '0'^ 



0.9 

o.e 



- r^ 



2 4 « BjQ., 2 

wd IN POUNDS (- 



Fig. 31 — Variation of Co with wd for cables of smootli exterior. 

The results of such a computation are shown in Table II for cable 
No. 1. For V > 1.5 knots the experimentally determined H is 64.0 de- 
gree-knots. The corresponding computed values of H, ranging from 67.4 
to 70.0 degree-knots, compare favorably ^ith this experimental value. 
Over the entire range of V, from 0.25 to 10.00 knots, the variation in // 
is about 4 per cent. This small vaiiation makes the assumption of Section 
3.2 that H is a constant for all V appear reasonable, especially since we 
can only hope to use this computation in a preliminary design before 
the hydrodynamic properties of a cable are established by experiment. 



Table II — Computed Values of A''^ , 
Cd , and H for Cable No. 1 



V (knots) 


^R 


Cd 


n (deg.-knots) 


0.25 
0.50 

1.50 

3.00 

10.00 


1.25 X 10^ 
2.50 
5.20 
5.80 
6.05 


0.935 
0.922 
0.965 
0.985 

1.000 


67.5 
70.0 
67.8 
67.5 
67.4 
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In Table III we have indicated computed high-velocity H values as 
a function of the unit weight w and the diameter rf of a cable. Table IV 
shows the computed high-velocity Co and H for various gauge piano 
wire. The American Steel and Wire gauge scale is used in this tabulation. 

Table III — Computed Values of the Hydrodynamic 
Constant H in Degree-Knots for Smooth Cable 





















Submerged 




















Weight in 
lb/ft 




















0.5 


0.75 


1.00 


1.25 


1.50 


1.75 


2.0 


2.5 


3.0 


4.0 


0.1 


54.5 


44.2 


37.9 


33.7 


30.6 


28.1 


26.2 


23.2 


21.0 


17.9 


0.2 


75.9 


61.1 


52.3 


46.4 


41.9 


38.5 


35.8 


31.6 


28.5 


24.4 


0.3 


91.7 


73.6 


62.9 


55.6 


50.2 


46.1 


42.8 


38.0 


34.4 


29.7 


0.4 


104.7 


83.9 


71.5 


63.1 


57.1 


52.5 


48.9 


43.5 


39.6 


34.3 


0.5 


115.9 


92.6 


78.9 


69.8 


63.3 


58.3 


54.4 


48.5 


44.3 


38.3 


0.6 


125.8 


100.4 


85.6 


75.9 


68.9 


63.6 


59.5 


53.2 


48.5 


42.0 


0.7 




107.6 


91.9 


81.6 


74.2 


68.7 


64.2 


57.4 


52.4 


45.3 


0.8 




114.2 


97.8 


87.0 


79.3 


73.4 


68.6 


61.3 


55.9 


48.4 


0.9 




120.6 


103.4 


92.1 


84.1 


77.8 


72.7 


65.0 


59.3 


51.3 


1.0 




126.5 


108.7 


97.1 


88.6 


82.0 


76.6 


68.5 


62.5 


54.1 


2.0 






153.3 


137.0 


125.0 


115.7 


108.2 


96.7 


88.2 


76.3 


3.0 








167.6 


152.9 


141.5 


132.3 


118.2 


107.9 


93.3 



Table IV — Computed Cd and // Values op Piano Wire 



Gftuge 
(Am. Steel & Wire) 




5 

10 
]5 
20 
25 
30 
35 
40 



Dia. (inches) 



0.009 
0.014 
0.024 
0.035 
0.045 
0,059 
0.080 
0.106 
0.138 



2.49 
1.91 
1.56 
1.39 
1.31 
1.24 
1.14 
1.02 
0.970 



n (deg- knots) 



10.7 
15.2 
22.1 
28.2 
33.0 
38.7 
47.2 
57.1 
67.1 



APPENDIX C 

Some Approximate Solutions for Laying and Recovery 

c.i Laying 

We assume that the tangential drag and the centrifugal forces are 
negligible. Then, since for laying ^ ^ g ir, (18a) by virtue of (21) be- 
comes 



(!•-) 



<^^ , A -3/. 

1- A sin 6 — cos 

as 



= 0. 



(61) 
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Let the origin of an x, y t^oordinate system be at the cable touchdown 
point (Fig. 8). Further, let x be the x coordinate of a point along the 
cable configuration and s the corresponding distance along the cable 
from the oiigin. If we define 

A = s ~ X, (62) 

then 

rfA ds dx , e ,^„, 

dy d\j dy 2' 

and 

^ = sin d. (64) 

ds 

By means of (63) and (64), (61) transforms to 

,-, _. dA rf"A 1 (dA)' (rfA)- 1 , , 

dy dy- 4 (dy) (dy) 4 

where we have in addition introduced the non-dimensional variables 

To = To/wh, 
A = A/h, 

y = y/h. 

Here h is the ocean depth at the touchdown point. Using the condition 
that 6 = at 1/ = 0, which implies dA/dy = at ?/ = 0, we get upon 
integrating (65) 

rfA ^ ^^^^ a I 1 - [%/{% + yW \" 



where 



2 — sin^ a ,„-^ 

7 - -^— . (67) 



The usual range of the critical angle a is between 10 and 30 degrees. Also 

0<-^^<l 
Therefore, we approximate the denominator of equation (07) by unity. 
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With the boundary condition that A = s - x = at y = 0, we thus 
obtain 

A(l) ^ ^ - X = tan ^ [' (1 - {To/{To + k)]')^ d^, (68) 

where S and X are the dimensionless values of s and x at the ship. 

Next we let 

to — s + a;, 03 = oj/k. 
Then we have 

dw .. I dL 
dy / dy' 

and, as can be seen from (06) through (68), 

^(1) =S-\-X= ctn^ f (1 - [fo/(To + $)]')"* d^. (69) 

J Jo 

For convenience we define u and R by 

To 



u = ^ 



n = 



n + r 

To 



1 + To" 

In terms of u and /e (68) and (69) become 

A(l) = % tan " C ^^—^ du, (70) 

Further, integration by parts gives 

Combining the above three equations and making the approximation 
(1 - H'')^ f^ 1, we find 



(l - l) A(l) + I tan= ^ ^(1) = (1 + To) tan ^ . 



(72) 



Thus A(l) and to(l) are related, and we need evaluate only one of the 
quantities numerically by means of equation (70) or (71) in order to 
compute both A(l) and w(l), and hence >S and X. 



DYNAMICS AND KINEMATICS OP SUBMARINE CABLE 1183 

The singularity at w = 1 makes the numerical evaluation of the inte- 
gral in (71) cumbersome. Therefore we consider the evaluation of A(l). 
But for convenience of numerical calculation we write (70) as 

A(l) = -rotan|£ (1 ~ u'Y^dfi^Y 

and integrate by parts to get 

Ad) = tan I ((1 - B^y - 1 % /; i^« ^^, ,u) . 

We note again that (1 — R'^y ^ 1. Further, essentially all of contribu- 
tion to the integral in this equation occurs near u = 1 because of the 
large value of 7. On the other hand, the values of 7^ which are of in- 
terest will normally be smaller than unity. Hence R, the lower limit, will 
normally be less than one-half, and thus will be outside of the region of 
significant contribution to the integral. Therefore, we can take the in- 
tegral to be a constant for a given a. Denoting this integral by n and 
combining these considerations we obtain 

A(l) = tan^ - ^ntan^n. (73) 

Finally solving for 8 and X from (68) and (69), we find 

sm a 

X = -J— + XTo, 

tan a 

which are a dimensionless form of (24) . For brevity wo have wi'itten 



1/2 

""six 



x=i? 



1 _ 2 ( ^ _ 1 K Ltn ^ + ^ " tan ^ 



2\7L 2V2 /J 22 2 

Since the integral n and the constant 7 depend only on a, the constants 
K and X are also functions of a only. We have evaluated n by numerical 
integration and have plotted the resulting values of k and X — k in Fig. 10. 

C.2 Recovery 

In the conventional recovery situation we have as boundary conditions 
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6 = at y = 0, 

(74) 
8 — —Ga at y = h, 

where a^ is the incidence angle of the cable at the ship (Fig. 8). With 
these boundary conditions, the development leading to (66) yields (26) 
for the relationship between shipboard tension 2\ and the incidence 
angle «« . 

To evaluate S and X, the values of s and x at the ship, we use (19), 
(20a) and (20b). Again we simplify by assuming Dt = PcVc^ =^ 0. Fur- 
ther, since for recovery —ir ^9^0, we may write | sin f | = —sin ?. 
This gives 



X 



Jo L^^os ^ + A sin'^ ^ 
Jo U 



exp 



f 

Jo 



sm T] 



cos 7j + A sin^ tj 



exp 



i 



sm t} 



dr, I d^. 



dv I d^, 



(75) 



_cos ^ + A sin^ ^ ^ Jo cos tj + A sin^ t} 

The dimensionless bottom tension To is computed from (26). The inte- 
grals appearing in (75) have been evaluated numerically. The results 
are shown in Figs. 14 and 15. 

APPENDIX D 

Analysis of the Effect of Ship Motion 

D.i Formulation of the Differential Equations 

To analyze the effect of ship motion on cable tensions, we use the model 
shown in Fig. 32. We assume the cable is a perfectly flexible and elastic 
string whose motion is planar. The distance L along the cable from the 
ship to the point of entry into the water is taken as constant, and the 
longitudinal damping as negligible. 



J_L 



1 i I I 1 P-* — I 



>J^ U 



AfR 



Fig. 32 — Model ii.sed for the analysis of ship motion tensions. 
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Unlike the solution of the basic stationary model, the complete solu- 
tion of this model is not simple. To make the problem tractable, we shall 
make further simplifying assumptions. Although these assumptions may 
seem reasonable, they must be ultimately justified by comparison of 
experience with predicted results. 

Force diagrams of a differential element of cable are shown in Figs. 
33(a) and (b) for the two regions, air and water respectively. The no- 
tation is 

p = longitudinal displacement of a point of the cable, 

q = transverse displacement of a point of the cable (in air), 

rj = transverse displacement of a point of the cable (in water), 

e = the stationary angle, i.e. the angle the cable configuration makes 

with the ship velocity in the absence of ship motion, 
tp = deviation from the stationary angle - positive in the clockwise 

direction, 
s = distance along the stretched cable, 
X = distance along the unstretched cable (in air), 
f = distance along the unstretched cable (in water), 
Wa = weight per unit length of cable in air. 




It^T h--^---^^ 




Fig. 33 — Diagram of forces acting ou a cable element in air and ia water. 
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Suimning forces along the directions t (tangential) and n (normal) showTi 
in Fig. 33, dividing by ^x (air) or Af (water) and letting Aa; — *■ and 
Af — > 0, we obtain the foUottdng equations of equilibrium* 



Air: 

Tipx — Wa cos {B — tp) = pciqtt cos (p — pn sin tp), 
T^ + Wa. sin {d — If) = pciqtt sin ip + pu cos (p), 
Water : 

Ttp^ + OjvSr ~ '^ cos (d — <p) = Pwiqii cos ^ — pu siu ip), 
T^ -]r w sin {6 — ip) = PciQu sin ^ + g^j cos ^). 



(76) 



(77) 



Here, pc denotes the mass per unit length of the cable in air. As is known 
from hydrodynamic theory, in order to accelerate a body through a 
fluid, one must change not only the momentum of the body but that of 
some of the surrounding fluid as well. Thus the body has a virtual or ap- 
parent mass in addition to its intrinsic mass. In the first (77), the equa- 
tion of equihbrium in the normal direction in water, we accordingly use 
p„, , given by 

Pvi — pc -r -^ a p 

as the intrinsic plus virtual mass per unit length of cable moving through 
water. The quantities d and p are the outer diameter of the cable and 
mass density of the water, respectively, and the quantity (7r/4)rf p is 
the virtual mass of a unit length circular cylinder movmg transversely 
through water. 

We take for the normal drag force per imit length 

D^ ^^V^\V^\. (78) 

Here Vn is the normal component of velocity of the water relative to 
the cable, i.e., 

Vn = V sin (d — <p) -\- lit sin tp — i^i cos (p, (79) 

and CDpd/2 is a constant. 

The quantities s and p are given by the following geometric relations 
which can be obtained from Fig. 33: 



* We use the subscript notation for differentiation throughout this section, 
e.g., fx =d'p/dx. 



DYNAMICS AND KINEMATICS OF SUBMARINE CABLE 1187 

s, = [(1 + Pif + V!;V (80) 

tanp = 7;r/(l + vt), (81) 

with similar expressions obtaining for the cable in air. 

The tension T in turn is given by the Hooke's Law or stress-strain re- 
lation 

T = EA{[{1 + p.f + q.'f - 11, (air) 

, 82) 

T = EA{[{l-\- vd + mY - 1}- (water) 

As we indicated in Section 4.1, we shall assume that the extensile rigidi- 
ties EA corresponding to complete restraint and no restraint to tmsting 
will give the limiting values of the ship motion tension. 

Equations (76) through (82) form a complete system in terms of the 
independent variables x or ^ and t. Formulating boundary conditions in 
terms of the coordinate x (or ^) is awkward. This coordinate is measured 
along the unstretched cable so that a disturbance applied at the ship is 
apphed at different x's as the cable is paid out. At the same time, if the 
velocity of the pay-out is small compared to the significant wave velocity 
of the cable then we can plausibly neglect the paying out effect. As will 
be shown subsequently, in the problem at hand there are two significant 
wave velocities, roughly corresponding to transverse and longitudinal 
motion. The first of these is of the order of 200 ft/sec, while the second 
is of the order of 5,000-10,000 ft/sec. On the other hand, the pay-out 
velocity is of the order of 10 ft/sec. Hence, we take the pay-out velocity 
to be zero. This allows us to use (76) through (82) without further trans- 
formations and to identify x and f as coordinates fixed in the translating 
reference frame. 

D.2 Perturbation Equations 

We assume that the motion is a small perturbation about the undis- 
turbed configuration of our model. To determine which terms of the 
differential equations are important in this ease, we adopt the following 
procedure. Let 

M - max[(Po + Pif + QoY, 

where Po and Qo are displacements of the cable at the ship, and Pi is the 
variation of the pay-out displacement from the mean. The quantity 
e ^ M /L will normally be less than unity, and for no ship motion will be 
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zero. We write 

Po-\-I\ = am, 

Qo - bgit), 

where f{t), g{t) are some bounded functions of time and a and h are 
constants. We assume that for given f{t) and g{t), T, p, and q vary 
analytically with e, namely 

r = To + eTi + e'T2 + ■ ■ ■ , 

q = egi + eq^ + ■ ■ ■ , (83) 

V = Po + epi + eVi + ■ ■ ■ , 

with counterparts for the submerged cable. The stationary transverse 
deflection is further assumed zero, and therefore the series for q contains 
no go term. Substituting, for example, (83) into (82) for air and equating 
like powers of e, we find 

T, = EApo. , (a) 

Ti = EApi^ , (b) ,g^^ 

T^=EA \V2. + T-r^— 1 ■ (c') 

Equation (84a) of this sequence shows that only longitudinal displace- 
ments are associated with stationary tensions, while (84b) indicates that 
for small ship motions cable tensions are independent of the transverse 
component of ship motion. To compute the effect of transverse motion, 
(84c) shows that terms of the order e in p and e in g must be considered. 
We assume further that 1 + po^ PtJ 1, since p^x is the order of magnitude 
of a strain. 

Equations (83) and (84) substituted into (76) yield with this approxi- 
mation 

Wa cos « — 0, (a) 



EA pdxx — pa Pott + Wa slu a = 0, (b) 

qux ^ qitt = 0, (a) 



(85) 



C2 

Pi.. --Put = 0, (b) (86) 

P2.T — — , Pitt = — 2 qutqu^ — quqixx , (c) 

Ci Ci 
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where 

C2 = Ta/pa . 

For non-zero Wa and a 9^ ir/2, (85a) cannot be satisfied. This is a conse- 
quence of the assumption of 50 ^ 0. With poi( = 0, equation (b) implies 
in turn that To = constant, which agrees with our model. For the sub- 
merged part of the cable, the equations do not yield a constant To and 
thus contradict the assumed model. However, on the assumption that 
the transverse motion is confined to a region near the surface, we con- 
sider To to be constant in the submerged part of the cable as well. We 
thus arrive at 

ijiff — 8i)i^ — yr]it — — rjitf = (a) 

pirf --Put = 0, (b) (87) 

Ci 



where 



C2 — To/pu, , 



, Cnp dV' 

8 = — ^ — - cos a sm a, 

CopdV" . 

as the differential equations governing the motion of the submerged 
cable. The constant Ci is the velocity of propagation of a longitudinal 
wave in the cable, while the constants C2 and C2 represent the propagation 
velocities of a transverse wave in air and water respectively. 

D.3 Solution of the Perturbation Equations 
We write 

p(o, t) = P,(t) + Pi(0, 

9(0, t) - Qoit), 
and take as boundary conditions 



1190 THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1957 

pi(0, t) = , (a) 

a 

vM t) = 0, (b) ^^^^ 

ei(0, t) - Qo/e. (c) 

That is, we apportion all of the longitudinal boundary motion to pi , 
and all of the transverse boundary motion to qi . Equations [(84b), (86b) 
and (87b)] then give the complete tension due to the longitudinal com- 
ponent of ship motion to first order. As mentioned in the text, this ten- 
sion is easily obtained from standard references, and is also the greater 
part of the ship motion tension. 

To determine the tensions due to transverse ship motion, we solve 
(86c) and (87c) for boundary conditions (88b) and (88c). In addition, 
we have the transition conditions 

Qi(L,0 = 7?i(0, 0, (a) 

9:. {L, t) = ^ir(0, 0, (b) (89) 

P2(x = L,t) = vM = 0,0, {^) 

p,AL,i) = pAO,t), (d) 

which follow if we assume that at the point of entry into the water the 
cable is continuous and the tensions are finite and continuous. 

We consider only the problem of the tensions associated with a har- 
monic steady-state transverse disturbance. Equations (86a) and (87a) 
show the transverse response to this disturbance to be independent of 
the longitudinal motion to first order. The first-order transverse motion 
in turn can be thought of as a forcing action on the second order longi- 
tudinal motion, as (86c) and (87c) indicate. This suggests the program 
we follow to compute tensions. Namely, we first determine the first- 
order steady-state transverse response, then the second-order steady- 
state longitudinal response which is excited by the first-order transverse 
oscillation, and finally, by (84c) the resulting tension caused by trans- 
verse motion. 

D.4 Transverse Response 

At the ship we assume a harmonic forcing function 

Qo{t) - A cos wi, (90) 

and we introduce complex exponential representation 
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qi = Re Qiix) e'"', 

m = Re ^i{f) e'"', 

where the factor e'"' will be henceforth suppressed. 

The solution of {80a) and (87a) for the steady state may then be 
written 

Qiix) = Bi exp — + B2 exp ( - — ), 

Co \ C2 / 

Fliit) = Fi exp (gif) + F2 exp (g^f), 

where the B's and F'a are complex constants and qi and 92 are the roots 
of the quadratic 

2 

q' — 5q — iciy + ^r, ^ 0. 
C2" 

Throwing away the root of this equation which corresponds to the in- 
coming wave in watei', we get 

H,{t) = F exp (qj^). 

where qi is the root corresponding to the outgoing wave. The three com- 
plex constants B\ , B2 , and F can now be determined from (89a), (89b) 
and (90) 

Si + 52 = A/e, 



— B\ exp — — Hi exp I — — - 1 



- q^F = 0. 



We note that Bi , Bi and F are proportional to the amplitude A of the 
forcing motion. 

D.5 Second-Order Longitudinal Response 

From the preceding results, the right-hand sides of the equations of 
longitudinal motion (86c) and (87c) can be computed. This computation 
for (86c) results in 
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K. 2a).r , •2to.v\ ^ , 

n sill — + Ti cos — 1 cos iwt 
C2 Ca / 

, / 'Zwx . 2a)a;\ . ., , 

+ r-i cos — Ti sm 1 sin lut 

\ Ca Ca / 

. 2t»)a; , 2aj.T; 
+ 7-5 sm — - + 7-2 cos , 

C2 Ch J 

where the r's, which are proportional to the square of the amplitude A, 
are 



(92) 



n 


= Re (S/ + B2-). 


Ti 


= Im (Bi' - B2'), 


n 


= Re (Bi' - Bi), 


n 


= Im (Br + B,'), 


?'b 


= 2Re B1B2 , 




2 2 
, era 



re -- r5 + 



A= 



[The quantity re will be used subsequently.] Similarly, for the right-hand 
side of (87c) we get 



VH 



/ _. ,j,„ - jjijf j = e '"^[{ai cos 2fff + a2 sin 2(rf) cos 2aii ,ggv 

+ (tti sin 2fff — a2 cos 2o-f) sin 2wi + aj, 



where 



Qi = — (p + ^'o-), 



ai = 



02 — 



Ql 



?1 



[( 



^1 COS (2/+^) + igi|^cos(2/+33)|, 



^ sin (2/ + 3) + I 51 I sin (2/ + 3^) 



]■ 



az=\F\ 1^^ J + I 91 
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and 



/ = arg F, 

g = arg (-51), 



g / ^ 27r, 



0^.^^. 



It is seen that expression (92) and (93) have terms of the form 

sin 2ut 



Fix) 



cos 2(jit 



(94) 



in addition to functions of x (or f) alone. In accordance with the idea 
that the first order transverse motion is a forcing action on the second 
order longitudinal motion, we take as solutions of (86c) and (87e) func- 
tions of the form 



Gix) 



sin 2hit 
cos 2(j}t 



to correspond to terms of the type given by (94) and functions of x (or f ) 
alone to correspond to forcing terms which are independent of time. This 
again gives linear differential equations which can be readily solved. For 
example, corresponding to the first term in (92) multiplying cos 2cct we 
have the assumed solution 



G{x) cos 2ut, 



and the differential equation 

/G, w^^ir/c,^^ 

dx' "^ ci^ 2 L\ci 

This has the solution 



G = L4.1 cos — 

L ci 



■f- A2 sin 1- 



8C2\ 



. 2a).'C , 2(j3X 

n sm ■ — ■ + r2 cos — 



. 2ioX , 2aiX 

n sin h Vi cos 

Ca C2 



]■ 



where Ai and A2 are undetermined constants. 

In this manner, the solution for the longitudinal motion can be ob- 
tained in terms of a set of con.stants. These in turn can be evaluated by 
means of the boundary and transition conditions on p2 ■ This evaluation, 
although straightforward, is very tedious. We shall omit the details of 
it here. The final result is 
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Air: 



e J 2 — 



4ciCi 



{[.,sin?^ + (. 



20)3: , / , ci . 2a>L 
sin + I rg ^ u sin 



Ci 2toL\ 2a)3; , Ci 1 o , 1 r ■ ^^^ 
-n cos cos + - ^6 cos 2iat + 7-3 sin — - 



Ci 2wL , ci . 2a)L\ 2a)a; , ci i ■ o » 
-ir- -\--ri cos — — + - re sin cos + -n} sin 2iot 



4 sii 



., , , 2a)a: 2ojL\ / . 2ux . 2uL\ | F 

+ -|r,(cos— -cos— l-r,(^sm— -sm-^j- — ^^^ 



[r.(c 



C3 



Waier: 



2 _ e'a)'£?A 



4CiC2 



2ti)L , . 2a)L 

— Ti COS + ra sin — - 

Ci Ci 



, ci . 2wL / . 2wL 2wL 

+ — sm n sm — re cos — 

C2 Ci \ Ci Ci /J 



sin 2w 



{-0 



+ 



. 2ajL , 2uL , ci . 2ajL / . 2l^L 

r2 sm + n cos + — sm ■ I re sin 

Ci Ci C2 Ci \ Ci 



+ ri COS 



2a)L 



Cl /J 



Ci 



cos2a,(i-ill- iT^r^e'^"^ 



Cl 



D.Q Numerical Results 

Since the r's are each proportional to the square of the amphtude A, 
the above results indicate that the transverse motion tension varies as A 
squared also. It is additionally a function of the frequency of ship motion 
oj, the forward mean ship velocity F, and the stationary tension To . The 
computation of the transverse ship motion tension for the laying situa- 
tion was carried out for cable No. 2. The results are shown in Fig. 34. 
Here we have denoted the transverse motion tension by Tg and have 
plotted Tg/A'^ against the period of ship motion t. Rather than the 
stationary tension 7\ , we have used the depth h, which during laying is 
du-ectly related to Tohyh= To/w. Fig. 34(a) is a plot of TJA'^ versus 
the period t - 27r/w for h = ^,2 and 3 nautical miles and for V = 6 
knots. Figure 34(b) is a plot of TJA^ versus t for 7 - 3, 6, and 9 knots 
and h = one nautical mile. 

For representative laying, for example at 6 knots with a ship period 
of 6 seconds into a depth of one nautical mile, Fig. 34 gives 
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TJA^ = 0.50 Ib/ft^ (twist unrestrained), 

Tg/A^ = 0.93 Ib/ft^ (tw-ist restrained). 

For an extreme value of A = 20 feet, we get therefore that Tq is be- 
tween 200 and 370 pounds. 

Additionally, by means of the above analysis, one can compute the 
rate of damping of a transverse disturbance after it enters the water. 
The results of this computation are shown in Fig. 18 and are discussed 
in Section 4.1. 
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Fig, 34 — Variation of the transverse ship motion tension of cable No. 2 
with the period of ship motion. 

APPENDIX E 

Tension Rise with Time for Suspended Cable 

E.i Formulation of the Solution of the Problem 

Let be the lowest point of the cable at time t after the suspension 
has begun (Fig. 35). We make the following definitions: 

h =^ depth at onset of the suspension, 

^1 = cable length from A to 0, 

tSo ^ cable length from ship at B to 0, 

.Y] = horizontal distance from A to 0, 
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Xa = horizontal distance from B to 0, 

5 = vertical distance from A to 0, 

To = cable tension at 0. 

If the cable is being paid out mth a slack «, then conservation of the 
total cable length gives the equation 



S, + S2 = 



k 



sni a 



+ (1 + t)Vt + cable stretching. 



(95) 




Fig. 35 — Coordinates for the analysis of tension rise when a cable is com- 
pletely suspended. 

It is assumed that there is no cable pulled from the bottom. The cable 
stretching we evaluate as in the example of Section 3.6, viz., 



cable stretching = (Si + S2) 



EA' 



This makes (95) read 



em a HjA 



(96) 



To obtain further relations for the unkno%vns appearing in (96), we 
assume that from the ship to point the cable configuration is a station- 
ary one governed by the equations developed in Section 3.G, while from 
points to A we assume that the cable configuration is a static catenary. 
These assumptions yield the following relations: 



S-i. = — sinh 0-, 
w 

62 = — : + K — , 

sm a w 



(a) 
(b) 
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X. ^i±i + x^, (c) 

tan a w 

5 = Tl (cosh a - 1), id) ^^^^ 

w 

X,-\-X, = -^ + Vt, (e) 

tan a 



"= T. 



(f) 



T. = To^- io{h + 5). (g) 



Here k and X are constants, defined and plotted in Section 3.0, which de- 
pend only on a, the critical angle corresponding to V. Equations (96) 
and (97) form a complete set of equations in the unknowns Xi , X2 , 
Si , Si , To , Ts , 8, and a. They can be reduced to a set which contain 
only the unknowns a and Tb : 



ipiia) sin a — t(p2{o) sin a — h[ 1 -\- — j-^t sin a ) = 0, (a) 

ip2{cr) 



^ii^) I „:. 



where 



T. = l + '-i!i^, (b) 



I = wh/EA, 
t = Vt/h, 
T, = T,/wh, 



(98) 



and 



<Pi = sinli a ~ a -\- (cosh a- — l) tan - — (X — k). 
cosh (T — 1 , . , 

<P2 = 1 + X + (T, 

tan a 

cosh ff — 1 , . , 

ws = r + snih o". 

sin a 

A graphical iteration method will be used to solve (98). First we solve 

^i(ff) sin a — tipi{a) sin a — S = (99) 
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for u by means of a nomograph to be described later. Next the quantities 

cosh <T fP'i{cd ■ 

are plotted as functions of o- for various a. These plots can then be used 
as follows to solve (98) for a given t. Solve (99) to obtain a^ . From the 
plot of <f>i{(r)/<p«{a) sin a compute 






Using the value h* for h, compute <ri from (99). With this value of tn , 
compute h* from 






etc. In this way a convergent sequence tro ,jti , ■ ■ - , (r„ is generated. 
Finally, from the plot of cosh o-/^2(o") obtain 1\ . 

The above iteration procedure sounds tedious. Actually, in most cases 
the iteration is not necessary because <r remains essentially independent 
of time. Thus, the solution of (99) by means of the accompanying nomo- 
graph will usually give the complete solution of the problem. 

E.2 Nomograph (Alignment Chart) for the Solution of Equation (99) 
The relations 



, lOipiic) d sin a 

xi = 0, X2 = d, x-i = 



1 + 10v>2(o") sin a ' 

■J. e <f>].W) sin a 

■iji ^ h, y-i = TKi y-i = 



(100) 



10 ' -^ 1 + \0ip2{<j) sin a ' 
where d is an arbitrary constant define parametrically three curves 

Vi = Vii^i), i = 1, 2, 3, 

which we imagine plotted on a cartesian (x, y) coordinate system. A set 
of values H, e, and o- determine three points {xi , iji) {i = 1, 2, 3) which 
he on these curves. If these points lie on a straight hne, it can be shown 
that they satisfy (99). 

On the left-hand sides of Fig. 36 we have plotted the curves given by 
(100) for various values of the critical angle a. The values of the parame- 
ters h = wh/EA and e, which describe the curves iji = yiixi) and 
y^ = ^2(3:2) respectively, are plotted on the indicated scales. Rather than 
indicate the values of a along the curve -1/3 = ^3(3:3) we have for con- 



DYNAMICS AND KINEMATICS OF SUBMARINE CABLE 



1199 




Fig. 36 — Nomograph for the solution of equation (99). 

venience made an auxiliary plot on the right-hand sides of Fig. :^G of 
7/3(0-) versus a, hut with numerical values of the ordinate omitted. 

In addition, we have plotted in Figs. 37 and 38 the functions 
\ftp2 (cosh a) and (^a/s^'a) -^i" « for various a. 

E.3 Numerical Example 



To illustrate the method of obtaining the tension rise mth time 
described above, we consider a numerical example for cable No. 2. The 
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Fig. 37 — Variation of 



cosb 



with 



values we assume for the parameters which enter the calculation are the 
following : 

e = 0.02, 

7 = 6 knots, (a ^ 12°), 

h = 6,000 ft, 

S = 1.2 X 10" lbs, 

h = 3.1 X 10^'. 

To solve (99), we connect on Fig. 36 the points e = 0.02 and H = 
wh/KA = 0.0031 with a straight edge and note the intersection with 
the intermediate 7/3 = ^/sC^a) curve for a = 12° (point A). We then locate 
the point on the i/3 versus o curve having the same ordinate (point B). 
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Finally, we obtain the root of (99), cr ^ 0.555 by reading off the cor- 
responding abscissa (point C). This value of a now serves as the starting 
point in the iteration procedure by which we find the tension correspond- 
ing to a given t. 

For example, for t = 1.0 (t = 600 seconds) we have the following se- 
quence of values 



a„ = 0.555 
0-1 - 0.580 
a2 - 0.580 



(>P3 sin a)/<pi S* = A [1 + (^ sin a) i/psl 

0.212 0.00379 

0.213 0.00380 



with a converging to the value 0.580. For a = 0.580, Fig. 37 gives 

cosh a/pa = 0.800. 
Hence, by (98b) for t = 1.0 

T, - 1.80, 
and 1\ = 1.80 wh or 7,600 pounds. 
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APPENDIX F 

The Three-Dimensional Stationary Model 

F.i Derivation of the Differential Equations 

Let "i, 5, k be unit vectors along the ^, v, f axes {Fig. 26) and7 a unit 
vector along the tangent to the cable configuration in the direction of 
positive s. As in the two-dimensional model, we take this to be opposite 
to the direction of travel of the cable elements along the configuration. 
With respect to the cable configuration the resultant velocity vector of 
the water is in the — z direction. We resolve this velocity into directions 
normal and tangential to the cable in the plane formed by i and t . The 
unit vector in the normal direction we denote by n, namely, 

^ = ~i+l^-?j . (101) 

I -1 + {i-t)t\ 

In analogy to the two-dimensional model we assume the normal and 
tangential drag forces depend only on the corresponding water velocity 
components. Thus, we take 

D, = ^a-nV)\ (102) 

Equilibrium of the forces acting on a cable element yields the equation 

T'^-\--i^-\-lDr-\-nD^ -]w = pM. (103) 

as as 

The vector a denotes the acceleration of an element of the cable as it 
moves at the constant pay-out velocity Vc along the cable configuration. 
It is easily shown that 

a =Vo'^. (104) 

as 

For convenience we introduce a second reference triad of orthogonal 
unit vectors?, u, and v as follows. The v vector is taken in the (|, t\ plane 
normal to7; the u vector is chosen equal to the vector product u X |. 
The angles ^ and 6 shown in Fig. 27 describe the orientation of the ( (, 
u, v) triad. In terms of these angles, we read from Fig. 27 the following 
table of direction cosines 





i 


J 


k 


7 


cos d COS ^ 


sin 9 


— cos 6 sin ^ 


u 


— sin d COS 1^ 


cos 8 


sin 6 sin i/* 


V 


sin \p 





cos i// 
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In the (t, u, v) system the vector n becomes for example 

u sm 6 cos ^ — j; sin ^ 
[siu^ 6 cos^ ^ + sin- \pY 

Imagine the origin of the {7, u,v) triad to traverse the cable at unit 
velocity. The triad during this traverse rotates like a rigid body with 
respect to the fixed (^, tj, f) frame. The rotation, which we denote by U, 
is seen from Fig. 27 to be 

ii = ] \p -\- V 9 = u cos d ^ -\- V d -^ 1 sin 9 4'- 

Here the dot denotes differentiation ^\ith respect to time, or since 
ds/dl = \ it may be interpreted as differentiation with respect to distance 
along the cable. The vector7 is a fixed vector of constant magnitude in 
the rotating (7, u,v) triad, hence 

^ = 5 X 7 = M0 - y cos e i//. (105) 

as 

From (101), (102), (104), and (105) we obtain for (103) 

{T - pcVc') Q^ - I cos H) +'i(^ ^ ^t) 

+ 9b^Z- (sin' e cos' \t + sin' ^)\u sin 6cos^|^ -v sin \^) (lOG) 

— w(u COS 6 -j- i sin 9) = 0, 

which gives the three scalar equations, (47). 
Further, let r(s) be the cable configuration, i.e., 

r{s) = Uis) +j v(s) +I-r(s), 

where $(s), ij(s), and f(s) are the f, t], f coordinates of a point s of the 
cable. Then 

7 ■: d^is) ,-zdvis) , 7 rff («) 
ds •' ds ds 

Forming the scalar product of 7 \vithi,j, k respectively, we get (48) of 
Section 7.1. 

In a 9, \l/, T space the solution trajectories of (47) are given by the 
solutions of 

3, do _ A(cos' ^ sin' 9 + sin" \p) cos rp sin 9 — cos 9 
(r-p.VJ^, ~ Dr/w - sin 9 ^~^' 

, _ 0. f/^ _ A (cos' }{/ sin' 6 + sin' xj/) ' sin ^ 
- P'^V'^)-^. cos 9{Dt/w - sin 9) ' 
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We see that the trajectories are periodic in both 8 and i^ with a period 
of 27r, and only a single region, say 

^ iA ^ 27r, 

need be considered. It is apparent that the straight hnes 

(1) ^ = 0, e^ a; (2) ^f- ^ 0, = a + x; 

(3) ^ = TT, 6 ^ 2ir - a; (4) ^^ = tt, 6 = ir - a; 

are solution trajectories which contain all values of T. Along other solu- 
tion trajectories in this region one easily verifies that 

T = ^"^"' 
w 



1 — exp / 



ho A[cos2 i/(d) sin^ 9 + sin^ i^ (6)]^ cos i/{d) sin 6 - cos 6/ ' 

where ^ = Tp(6) is obtained from the solution of 

dip _ [cos d — A(cos^ i/' sin^ 9 + sin^ tf) cos ^ sin 6] cos 6 
(16 A(cos^ ^ sin^ B + sin^ ^)' sin ^|/ 

From the definitions of \p and 6, it follows that the lines (3) and (4) are 
physically identical with Hnes (1) and (2), and represent straight-line 
laying and recovery respectively. Likewise, the expression for 7' shows 
that any non-straight line trajectory with zero bottom tension is bounded 
by PcVc/w- Hence, as in the case of the two-dimensional model, we con- 
clude that if the tension is somewhere greater than pcVc/w and the bot- 
tom tension is zero, the only possible stationary configuration is the 
straight line lying in the plane of the resultant ship velocity and gravity 
vectors, and making the critical angle a with the horizontal. 

F.2 Perturbation Solution for a Uniform Cross Current 

At the outset we assmne the tangential drag force to be zero. This 
gives by (49) 

T = w{h + r,), (107) 

where h is the total ocean depth. Furthermore, we take peVo to be zero. 

If the angle <p (Fig. 26) is small compared to unity, we assume that 6 

and ^ will vary only slightly from the values they would have if the 
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Upper, cross-current stratum extended all the way to the ocean bottom. 
That iy, we take 6 and ij/ to be of the form 

6 = a' -\- 8, 

(108) 

where a' is the stationary incidence angle corresponding to the velocity 
V, and 6 and jj/ are assumed small compared to unity. 

Substituting (48b), (107), and (108) into (47a, b) and retaining only 
linear terms in 6, tp and their derivatives, we get the linear first order 
equations 

(h + ,,) ^ + (2ctn-a' + 1)9 =0, (a) 

dv 

: (109) 

(h + Tj) ^ + cscV',^ - 0. (b) 

0.7) 

Because in the lower stratum the cable is a straight line parallel to the 
path of the ship, we have as boundary conditions: 

, ^ - h',{'- = " - "' (110) 

whore h' is the depth of the upper, cross-current stratum and a is the 
stationary incidence angle corresponding to the velocity V, 
The solution of (109) for the boundary conditions (110) is 

h - k'Y , 



k + n. 

(Ill) 



where 

fi = (2 ctn'a' + 1), 

p - cscV, (112) 

Aa = a — a'. 

Equation (48) for the space-coordinates ^, t], and f of the cable in turn 
can be written to terms of first order in the form 

-7- = etna — wsc a , 
dv 

d ^^^^^ 

-_ = — Ictna. 
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Substituting (111) into (113) and integrating under the condition that at 
^ = 0, ^ = and f = 0, we find 

^ , (k- h')Aa , , r 1 1 1 

^ = .etna + -^^-^ CSC . L^^qr^pr - ;T^J . ^^^^^ 



^ , ,{h-h') 
f = etna -. -T- ip 



1 1 



J- 



These equations describe the space cur\'e formed by the cable in the 
cross-current stratum. 

To determine the distances d and e (Fig. 32), we transform (114) for 
the cable configuration to coordinates $' and f' oriented along the ship's 
path and normal to it respectively by means of 

^' = I cos <p — t f^hi <p, 

f ' = I sin <p -\- ^ cos <p. 

The result to terms of the first order is 

. ' Ah- h')Aa 2 ,[ 1 1 1 

^ = .etna + ^-— ^csc a [^^ ^ ^y., - ;^J , 

/ , , Ctn a'jh - h') r 1 1 1\ 

Letting ij = —h' and denoting the corresponding values of J' and f' 
by — rf and —e respectively, we obtain (52). 
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